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ABSTRACT 


Major components of ices on interstellar grains in molecular clouds - water 
and carbon oxides - occur at various optical depths. This implies that selective 
desorption mechanisms are at work. An astrochemical model of a contracting 
low-mass molecular cloud core is presented. Ice was treated as consisting of 
the surface and three subsurface layers (i.e., sublayers). Photodesorption, reac¬ 
tive desorption, and indirect reactive desorption were investigated. The latter 
manifests itself through desorption from H-|-H reaction on grains. Desorption of 
shallow subsurface species was also included. Modeling results suggest the exis¬ 
tence of a ’’photon-dominated ice” during the early phases of core contraction. 
Subsurface ice is chemically processed by interstellar photons, which produces 
complex organic molecules (COMs). Desorption from the subsurface layer re¬ 
sults in high COM gas-phase abundances at Ay = 2.4-lOmag. This may 

contribute toward an explanation for COM observations in dark cores. It was 
found that photodesorption mostly governs the onset of ice accumulation onto 
grains. Reaction-specihc reactive desorption is efficient for small molecules that 
form via highly exothermic atom-addition reactions. Higher reactive desorption 
efficiency results in lower gas-phase abundances of COMs. Indirect reactive des¬ 
orption allows for closely reproducing the observed H 20 :C 0 :C 02 ratio toward a 
number of background stars. Presumably, this can be done by any mechanism 
whose efficiency hts with the sequence CO > CO 2 >> H 2 O. After the freeze-out 
has ended, the three sublayers represent chemically distinct parts of the mantle. 
The likely Ay threshold for the appearance of CO ice is 8-10.5mag . The lower 
value is supported by observations. 

Subject headings: astrochemistry - molecular processes - ISM: clouds - ISM; 
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Introduction 


Interstellar ices occur on the surface of dust grains in dark interstellar clouds. Cold and 
dense molecular cores form in the clouds via gr avitational contract ion. The cores evolve 


into protostars over a time-scale of order IMyr flLee &: MversI 1199911 . During this evolution, 


molecules settle onto the dust-grains. In dark cores, almost a complete freeze-out occurs. 
An ice mantle, probably hundreds of monolayers (ML) thick, acc umulates onto a typical 


grain. Surface reactions are a major route for molecule formation flPickles fc Williams 


197 


in interstellar ices. Along with accretion, this process dehnes the initial composition for all 
subsurface ice layers. Gas and surface chemistry are influenced by thermal and non-thermal 
desorption mechanisms. Selective desorption mechanisms, whose efficiency depends on 
surface species’ adsorption energy Ed, is of particular importance because of its ability to 
affect the relative proportions of abundances for solid species. Non-selective mechanisms 
also may have an effect on these proportions by determining the starting time and length of 
the ice formation epoch. These temporal parameters, in turn, affect the physical conditions 
at which the formation of ice occurs. 

The basic source of motivation for this article was the incomplete understanding 
on the role of desorption processes in dark clouds. One of the last mod eling studies 


on the relative signihcance of several desorption mechanisms was that of 


Roberts et ah 


(120071 ). Since then, several advances in astrochemical modeling have been incorporated 


in the models on a common basis. They in clude, hrst, real time-dependent cloud core 


models with changing physical parameters flBrown et ah 
and, seco nd, the consideration o 


1988 


Garrod fc Herbst 


2000h 


models) (IGuppen fc Herbst 


2007 


subsurface layers of interstellar ice mantles (three-phase 
). Additionally, reliable information from experiments 


or theoretic al studies has been published on desorption mechanisms, such a s reactive 


desorption flGarrod et ah 


20071) and photodesorption (e.g. 


Oberg et ah 


2009ah . 
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These new data 


( Garrod &: Pauly 


Chang &: Herbst 


l ave been implemented in 


201 


2014 


Garrod 


2018al: 


Du et ah 


lO recent astrochemical models 


2012 


Vasvunin fc Herbst 


2013aib 


e.g.). Reac 


dve desorpt i on has been analyzed in d etail as a non- 


selective mechanism (iGarrod et ah 


2006 


2OO7I 


Vasvunin fc Herbst 


other papers describe it as a selective mechanism fjPn et al 


2012 


2013b), only, a 


Rebonssin et al 


thoug h 


20i4 . 


All these authors have used a two-phase gas-grain model. In the case of desorption by 
interstellar and cosmic-ray-induced photons, a standard yield of 0.001 desorbed molecules 
per photon is typically used in the abovementioned papers. Instead, values closer to 
0.002 or 0.003 are supported by experiments (see section Finally, desorption by 
the exothermic H-l-H reactions on grai ns was identihed as p robably the most important 


desorption mechanism in dark cores by 


Roberts et al. 


(120071 1. yet, it has not been included 


into any of the recently published models altogether. 


Following the above discussion, the aims of this paper are to analyze (1) the signihcance 
and selectivity of higher-yield photodesorption (section [3^ : (2) the effects of selective 
reactive desorption (section : and (3) the importance of indirect reactive desorption, 
manifested by the H-l-H reaction fsection 13.411 . The analysis has been performed with a 
time-dependent three-phase model. This allows to offer a fresh view for aims (2) and (3), 
and is a requirement for an up-to-date investigation for (1). 


1.1. Ice physical model development 


In order to perform an up-to-date analysis of the desorption mechanisms, an advanced 
astrochemical model is required. The models employed for the physical description of 
interstellar ices have often been rather limited. Ices accumulate in mantles with a thickness 
of around 10^ monolayers (ML ). The (inert) bulk of the icy mantles was hrst considered 


by 


Hasegawa fc HerbstI fjl993a|) . Yet many studies continue to rely on two-phase (gas 
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and surface) models. Any reactions in ice mantl es have to be hampered by the obvious 


difficulties of molecule movement in ice lattice. 


Kalvans fc Shmeld ( 2010 . 2013 ) have 


investigated the possible result of mantle reactions that occur on the surface of cavities 
within the mantle. The cavities contain up to a few per cent of ice molecules at a given 
time, which mean s that mantle r eaction efficiency is at most, only tenths of that for 


surface reactions. 


Gar rod! ()2013al) employs a different approach ~ mantle molecules are 


more strongly bound into the ice, and thus move more slowly, with the reaction rates 
affected accordingly. These are, literally, three-phase models, i.e., the sub-surface mantle 
was considered as a fully mixed phase with no distinction between ice monolayers. This is 
a major approximation because it has been shown that ice composition is depth-dependent 


(iGarrod fc Paulv 


20111 ). 


An effective option for ice-mantle modeling is the Monte Garlo technique. This 


approach initially suffered from a limi 


computationally intense calculations fiGuppen &: Herbst 


ed set of surface specie s and reactions 


2007 : 


Guppen et ah 


compu tational cost places some restrictions for more recent models, too (IVasvnnin fc Herbst 


owin g to the 


20091). The 


2013af) . The approach has been used to validate the results of the macroscopic rate-equation 
models. One of the main conclusions is that rate-equation models may successfully 
reproduce chemistry on interstellar grains, only, if molecule binding energies Ef, are 
sufficiently high, e.g. Ef, > 0.5i?o, where Ejj is the adsorption (or desorption) energy of 
the particular species. High Et means that for reactive specie s accretion from the gas 


is faster than the time to scan the surface of the whole grain (IVasvnnin et ah 


20091). A 


study of Mo nte-Garlo random-walk 
published by 


Ghang fc Herbst 


(I2014J) . 


iechnique application for bulk-ice chemistry is has been 


In addition to the investigation of desorption mechanisms, a secondary aim is to 
develop further the ice-mantle description with the addition of a depth dimension for the 
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subsurface layers. An approach similar to that presented by 


GarrodI fl2013a() was us ed. The 


rate-co efScient calculation for mantle reactions is revisited because the approach of 


Garrod 


fl2013aj) has been taken directly from the description of surface chemistry. For surface, an 


underlying assumption is that molecules can freely migrate. Instead, for mantle chemistry 
the assumption used here is that molecules cannot move, they are mostly frozen in place. 


1.2. Significance of desorption of ice molecules 


Many surface chemistry studies have been performed without a proper reproduction 
of the observed abundances for important ice species. In particular, the production of 


GO 9 has been a major 


problem (iRnffle fc Herbst 


2 nnibl ) that has recently been solved by 


Garrod fc Paulvi (1201 ll) . The ‘appearance’ of GO 2 ice at a threshold Ay value was linked 


to the molecule binding energy, while residing on surface. The synthesis of GO 2 is also 
certainly connected to the rate of accumulation of oxygen and carbon in ice, which is largely 


governed by non-thermal desorption mechanisms (jLeger 


1983D. 


Ice composi 


Leger et al. 


1985 


i on as a function of desorption has been repe atedly re-visited (e.g.. 


Willacy fc Williams 


1993 


Roberts et al 


20071 ). Several major advances 


in the held require a renewed dedicated investigation. These are the rec ognition of the ro le 
of th e interstellar UV radiation held (ISRF) in regulating ice accretion fl Garrod fc Paulv 


2008 


2011 ) , the possibi 


Oberg et al. 


ity of s ubsurface mole cule photodesorptio n flAndersson &: van Dishoeck 


of reactive desorption flDu et al. 


2009bl). codesorption flOberg et al 


2012 


Vasvunin fc Herbst 


2009bl). an d the possible importance 


2013bh . 


Major interstellar ice components have been detected at diherent ISRF extinctions. 
The respective threshold Ay fAtv,) for H 2 O, CO 9 , and GO ices are app roximately 3.2, 4.3, 


and 6.8mag, respectively fjWhittet et al 


2001 


20071: 


Bergin et al. 


2 OO 5 I ). These values show 





































































a dependence on the desorption energies Er, for these species. Photoprocess, the principal 
desorption mechanism at low Ay, nsnally is assnmed to be non-selective. This is cansed by 


the non-thermal n atnre for this desor ption mechanism and codesorption of different species 


in an ice mixtnre flBertin et ah 


20131) . This disc repancy hampers the attempts to have a 


full ag reement between models and observations flGarrod fc Panlv 


2011 


Vasvnnin fc Herbst 


2013af) . 


Thus, a major aim of this paper is the investigation of reactive desorption and 
photodesorption that are Ay-dependent, and likely govern the accumulation of ice layers. 
These mechanisms have previously often been assumed being non-selective. The study 
involves empiric or semi-empiric equations for a proper calculation of desorption rate 
coefficients. 


1.3. Desorption mechanisms 


In the following paragraph, I shortly review desorption mechanisms, used in the present 
study. 


Evaporation from grains, which are in thermal equilibrium with the gas, is known to 
be an effective desorption path for molecules with very small Ejj (H, H 2 , He), or at higher 
temperatures. It is generally not sufficient for explaining gas-phase m olecular abundances in 


19831 ). After a protostar has 


dark clouds, which require other desorption mechanisms (iLegerl . 
been formed in a collapsing cloud core, gas and dust temperatures grow, and evaporation is 


among the dominant mechanisms flBrown et ah 


19881 ). 


Desorption hy the ISRF photons has often been omitted in models of dark cloud cores 
because the radiation held is very low at high extinctions {Ay > lOmag). However, the 
freeze-out stage for volatile species onto grains is largely governed by the ISRF. As the 
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dense cores evolve, the ice-forming molecules either disappear into the central body or 
are sublimated by the energy released by the newborn protostar. These processes have 


been fully reflected in astrochemical models only relatively recently (e.g., 


Garrod &: Herbst 


20061 ). One of the results is that the lifetimes of the ice mantles are largely determined by 


interstellar-photon-induced desorption. 


Cosmic ray particles that interact with the abundant molecular 


an internal UV radiation held within the clouds fjPrasad fc Tarafdar 


hydr ogen generate 


19831) . The resulting 


secondary photons are an important int ernal ioniza 


desorption of ice molecules from grains flDulev et al 


don source, and are able to induc e 


1989 


Hartquist &: Williams 


19901) 


Desorption by photons has traditio nally been treated a s a non-selective mechanism 


Roberts et al. 


20071) . The use of mole cule-specihc 


with uniform yields for all species (e.g., 
yields ha s become possible with the publication of reli able experimental results flOberg et al 


2009a 


u 


Rahr fc Baragiola 


2012 


Favolle et al. 


20131) . However, the phenomenon of 


codesorption may mean that this 


with similar p 


(IFavolle et al. 


lysic a 


2013 


Berlin et al. 


process is indeed non-selective, at least for species 


properties (Berlin et al. 


20131) . Wavelength-dependent studies 


20131) have revealed that the yield is different for 


conditions at cloud edges and cloud center, where ISRF and cosmic-ray-induced photons 
dominate, respectively. 

Whole-grain heating by cosmic-ray iron nuclei and subsequent ice sublimation 
has been recognized as one o f the principal desorption mechanisms in dark cloud cores 


fiLeger. Jura, fc Omont 


19851) . These auth ors consider evapora. t ion of C O from grains heated 


above 25K. In another classical paper by 


Hasegawa fc Herbst 


fll993al) . the evaporation is 
approximated as occurring on grains at T^arm = 70K for a time fwarm = 10“^s after each 
hit of iron cosmic-ray nuclei. The latter approach has often been used by subsequent 
astrochemical models. 
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Indirect reactive desorption by the excess energy of exothermic reactions on grains ha s 


been 


1972a 


bund to be a potentially important mechanism f Williami 


1968 


Watson fc Salpeter 


bl). It is driven by ionizing radiation helds that generate radicals and ions in gas and 


ice. The subsequent surface reactions release the energy, initially delivered by photons or 
cosmic rays. This type of desorption is tied to the radical content in the cloud, which is 
determined by ionizing radiation flux in it. Thus, reactive desorption is an ^dy-dependent 
process up to a certain depth, when cosmic-ray induced effects st art to dominate. This 


mechanism can be very important regarding the formation of ices flGarrod &: Paulv 


201 ih . 


As a separate 
was considered by 


Dulev & Williams 

(1993 

) and 

Willacv et al. 

(1994) 


(I 1994 J) . In addition to the 


ejection of the newly formed H 2 molecule, excess energy transferred to ice lattice may result 


in the d esorption of a near 


few eV flRoberts et al. 


12007 


Dy heavy atom or molecule. The excess energy is of order a 
, much less than that of UV photons. The yield is currently 


unknown, and is more constrained by its effect on modeled molecular abundances in clouds. 
This process has been assumed to be effective only for species with low Ed, e.g., for CO 
and lighter species, and thus it is a selective desorption mechanism. 

Other possible desorption causes include molecule ejection locally by cosmic rays. 
X-rays (important in protostellar disks), and chemical explosions. Of all these mechanisms, 
reactive desorption, desorption by the H-fH reaction and whole-grain heating are truly 
selective. The latter is rather ineffective. Experiments have sho wn that photodesor ption 


fISRF and cosmic-rav-induced) disph 

Favolle et al. 

2013; 

Berlin et al.l 

2013) 


ays some selectivity, too fjOberg et al. 


2009b 


201311 . Thus, direct and indirect reactive desorption in 


addition to photodesorption can be considered as selective and important during the 
formation of interstellar ices. These processes are investigated further in this paper. 


In section [2] a description of the chemical and physical model is provided. Each 
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desorption process has been studied individually in Sections 13.2113.31 and 13.41 respectively. 
In section 13.51 the results of a model that combines elements of previous simulations are 
investigated. Concluding remarks are summarized in section HI 


2. The model 

2.1. Sublayer approach in ice modeling 


A modihed code of the ALCHEMIC astrochemistry program ( Semenov et ah 2010 1 
was utilized to perform chemical kinetics calculations. Ice-related processes have been 
added (e.g., desorption mechanisms, the formation of subsurface ice mantles, and molecule 
diffusion in ice), and a real time-dependency, with changing temperature, density. Ay, and 
grain size, has been introduced into the model. 

An important feature of the model is the treatment of subsurface ice mantle. It has 


been recently recogni zed that mo. 


of thermal diffusion flOberg et ah 


ecules i n ice lattice may cha nge their location as a result 


2009d: 


Favolle et al 


2011a|l . Meanwhile, reactions are 


taking place directly in ice lattice, too, evidenced by the formati on of different molecu 


es 


in ices in the presence of radicals or ot 

ler reactive s 

pecies (e.g.. Gerakines et al. 1996; 

Oberg et al. 201 lb: Linnartz et al. 

o 

CM 

Noble et al. 

2013). This has been reflected in a 

recent astrochemical model by ( 

GarroJ 

2013a 

) that includes reactions among molecules 


in subsurface ice. More rigorous methods are those based on the Monte-Carlo techni que, 
although they rnay ha ve a high computational cost (e.g.. 


Chang fc Herbst 


201J). 


Vasvunin fc Herbst! 


2ni3at 


In order to include subsurface ice chemistry into the simulation, ice was described as 
consisting of four layers - the sur f ace an d three mantle layers. Such an approach lacks the 


resolution of 


Hasegawa fc Herbst 


fjl993bl) and similar models that consider ice monolayers. 
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Calculating chemical reactions in each monolayer for each time-step of the simulation is 
and enormous computational task. This is probably unnecessary, given that interstellar ices 
likely are characterized by irregularity. 


Ice accumulates onto irregular grains, it is amorphous, and may contain inhomogeneities, 
such as pores, large organic species, PAHs, or even small gra ins. Photodissociation 


products almost instantly 


Andersson fc van Dishoeck 


migra te across several monolayers flAndersson et ah 


2006 


20081) . A similar effect is possible for the products of exothermic 


binary reactions. In such an irregular ice, it is hard to define the boundaries and properties 
of a monolayer. Additionally, a radical in ice may react with any of its neighbors, some of 
which are in different monolayers, and it is unclear, which monolayer should the product 
species belong to. Most of these irregularities are of no concern, if ice layers of ten or more 
monolayers thick are considered as whole entities. 


Thus, it can be physically adequate to employ ice model that describes the ice as 
consisting of several subsurface layers instead of hundreds of monolayers. For short, I refer 
to these former layers as ‘sublayers’ hereafter. Each sublayer may contain a few up to 
dozens of monolayers. Molecules are assumed to be intimately mixed within a sublayer. 
Chemical reactions are permitted between same-sublayer species, only. Molecule exchange 
via diffusion can occur between adjacent sublayers. 


The number of sublayers can be chosen depending on the aims of the work. A minimum 
number of two is necessary to represent the two m ost distinct ice environments - pola r 


H^O- dominated and non-polar CO-dominated ices fiSandford et ah 


1988; 


Tielens et al 


19911 ). Three or more sublayers would be more adequate to perform calculations with 


the resolution required for a model that is able to fully account for ice observations 


in dense cores. This is much less than the hunc 


considered for some of the other bulk ice models flGarrod fc Paulv 


reds of monolayers that have to be 


2011 


Taquet et al. 


2012 
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Vasvunin &: Herbst 


2013al) . The low number of sublayers allows the explicit consideration 


of subsurface reactions in a depth-dependent manner. Three was chosen as the number 
of subsurface layers for the current research. They were numbered, beginning from the 
surface, with sublayer 1 being the outermost and sublayer 3 the innermost layer (connected 
to the grain itself). The outer surface is a separate, ‘zeroth’, layer. The formation of ice is 
described in section 12.41 


2.2. Physical model 

The rates for a physical phase-change or ice layer transition processes and chemical 
reactions are calculated via classical chemical kinetics equations. The rate coefficients for 
specihc processes are given in sections below. In cases, when phase-change occurs, we use 
the denominations /O and / to indicate the initial and hnal phases, respectively. In this 
regard, the ‘phases’ are gas {g), surface (5), mantle sublayer 1 (Ml), sublayer 2 (M2), or 
sublayer 3 (M3). 


2.2.1. Cloud conditions 


The model conside rs a gas clump in a cloud core in isothermal collapse. This approach 


has been adopted from 


Garrod fc Herbst 


was calculated according to 


Garrod fc Paulv 


(120061). T he temperature T, dependent on Ay, 


(1201 ih . The time-dependent results may also 


be interpreted as being dependent on the distance from the center of the contracting core. 
In such a case, earlier ti mes represent gas clu mps that reside further from the center. 


similarly to the model by 


Garrod et ah 


(l2008h 


The physical conditions were adopted for a 2Mq core and initial density uh = 
3 X 10^cm“^. The density increases to lO^crn”^ over a period of IMyr, in line with 
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Brown et alJ (119881) and iNeiad et al.l fjl990l) . It was assumed that a screen of diffuse gas with 
Ay = 1 surrounds the contracting core. The visual extinction to the center of the clump 
was calculated self-consistently, and grows from belo w 2 to over 200mag . The p hysical 


conditions in the cloud are similar to those figured in 


Vasvunin fc Herbst 


( 2013al ). The self¬ 


and mutual-shi 

elding 

from 

Lee et al. 

(19961 


data from 


ding of CO and II 2 has been included with the use of the tabulated data 


Li et al 


(Il996h. T he same was done for the N 2 molecule with the use of tabulated 


0201 31). 


The abundances of the chemical elements have been adopted from 


Garrod &: Herbst 


fl2006l) . Because a partially shielded molecular core is modeled, all elements, except for 


hydrogen, initially are in neutral atomic form in the gas phase. Hydrogen is divided between 
molecular H 2 (99.9%) and atomic H (0.1%). 

Hydrogen i onization rate by cosmic rays was assumed characteristic for dark clouds, 
1.3 X 10“^^s“^ QTomasko fc Spitzeiill968l) . The flux of cosmic-ray induced photons was 


taken to 


— 1 —2 
s cm 


pe Frrph = 4875s ^c m and the flux of heavy cosmic rays FpeCR = 2.06 x 10 


-3 


Roberts et al 


20071) . Finally, the assumed flux of ISRF far-ultravi olet pho 
photodesorption calculations was Fisrf = 1-7 x 10®s“^cm“^ at Ay = Omag flTielensI 


;ons 


200,511 . 


or 


2.2.S. Grain properties and grain-related effects 


The ’standard’ reaction network of 


Laas et al. 


(1201111 was used for gas-phase reactions. 


gas-grai n interactions and solid-phase rea ctions. This database, in turn, was based on the 


work by 


Garrod. Weaver, fc Herbst 


energy was taken frorn 
was taken 290K from 
added from 


Katz et al 


Roser et a. 


Vasvunin fc Herbst 


(I 2 OO 8 II . Three changes were introduced - H 2 adsorption 


1999 1. energy barrier for the GO -|- O surface reaction 


(I2001I1 . and the gas-phase reaction GH 3 O -f GH .3 was 


(j2013bl) . Following the experimental results of 


Katz et al. 
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fll999[) . quantum tunneling for diffusion and reactions was not considered. For reactions and 


physical processes, the model considers grains with a radius a = O.l/xm and an additional 
ice thickness b. The abundance of grains Ug = 1.3 x was calculated for spherical 

particles with a density of 3g cm“^ that contain 1% of cloud mass. It was assumed that 
the grains are rough, having Ng = 1.5 x 10® adsorption sites on the surface. This number 
was kept constant during cloud evolution under the assumption that surface roughness is 
smoothed out by the accumulating ice layers, while the size of the grain increases. Where 
necessary, it was assumed that a ‘cubic average’ ice molecule has a size of 5^ = 3.7 x 10“®cm. 
Maximum ice thickness in monolayers B = b/bm, achieved by this model, does not exceed 
180. 


2.3. Gas and gas-grain interactions 


Gas-phase processes include binary reactions, dissociation and ionization by interstellar 
and secondary photons, and ionization by cosmic rays. The grain albedo was taken to be 
0.5. Gas-grain interactions consist of (dissociative, non-adsorptive) recombination between 
positive ions and negatively charged grains, molecule accreti on (freeze-out), and d esorption. 


These processes were included as in the ALGHEMIG model fjSemenov et ah 


2010 h for grains 


with a radius a + b, with b calculated self-consistently. The sticking coefEcient was taken 


unity for heavy species and 0.33 for hydrogen atoms and molecules flBrown fc Gharnlev 


1990). 


In the reference, or ‘Standard’, model four desorption mechanisms, (/ = g, fO = S), 


have been included. They are evaporation at the e quilibrium grain 


evapor ation at 70K, induced by whole-grain heati ng. flLeger et ah 


1993a), p 


lotodesorption by i nterstellar photons flTurnei 


photons flRoberts et al 


1985 


t emperature T, 


Hasegawa fc HerbstI 


19981) . and by cosmic-ray induced 


20071) . Additionally, reactive desorption was considered, assuming 
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that 1.0% of surface reaction products are released into the gas fiGarrod et al 


200711 


These mechanisms for the Standard mode 


models by other autho rs fiGarrod fc Paulv 


Albertsson et ah 


2018h . 


have been c 


losen to be in line with recent 


2011 


Garrpd 


2013a; 


Vasvunin &: Herbst 


2013a; 


The desorption 


1 vield for the secondarv P 

lotons has been measurer 

or calculated 

(Andersson & van Dishoeck 

2008; 

Obere et al. 

2009a 

b), while the 


yield for most other species is unknown. In order to avoid unnecessary bias, an uniform 
yield of Fcrph = 10“^ desorbed molecules per photon was assumed for all species in the 


Standard model. There ar e experimental results t 


lat indicate a 


for higher-energy photons (IBahr fc Baraeiola 


2012 


Favolle et al. 


dependent on the absorption spectra of a particular species flBertin et al. 


righer desorption yield 


2013h. a 


houg h it is 


2013h . For 


consistency, it was assumed that Fcrph = fhsrf- The rate coefficient for photodesorption is 

7r(a &)^FphTph 


^pd — 


iV, 


( 1 ) 


where ‘ph’ denotes either ISRF or cosmic-ray-induced secondary photons and a + b 


is the time-dependent grain radius 1 section I2.2.2p . The hxed value of Ng means that 
photodesorption efficiency grows with the number of adsorbed molecules (up to 3ML, 
see b elow). This takes into a ccount the localized nature of photon absorption in a grain 


(e.g.. 

Dulev & Williams 

1988 

) and the dependence of photodesorption yields on ice depth 

(Andersson et al. 

2006; 

Andersson & van Dishoeck 2008: Obere; et al. 2009b: Bertin et al. 

201sT 

in a simple and natural way within the model. 


Photodesorption has been attributed to all sublayer species, given that the number of 
above monolayers does not exceed two. The desorption rate from a sublayer was divided by 
the number of monolayers in that sublayer. This accounts for the fact that only the upper 
ML has been exposed to desorption. Thus, desorption from a depth of up to three MLs 
is possible. This approach means that a number of ice monolayers can accumulate onto 
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grains at higher Ay values, earlier in the evolution of the cloud core, before photodesorption 
reaches its maximum effectiveness. 

In addition, desorption from the molecular hydrogen formation reaction on grains 
was included in the model for case investigation fsection l3.4p . The selectivity of ISRF, 
secondary-photon induced photodesorption, and reactive desorption was investigated in 
section |3l too. 


2.4. Physical description of ice 


S.4-J. Properties of ice layers 


It was assumed that ice consists of a total of four layers - the surface (abbreviated 
S), and three subsurface mantle layers, or sublayers 1, 2, and 3 (Ml, M2, and M3, 
respectively). Ml lies beneath the surface, while M3 is bound to the grain core. Molecules 
within each of the layers are completely mixed together. This approach is in line with 
recent similar models that consider subsu rface ice chemistry with the rate-equation method 


(iGarrod 


2ni3al: 


Kalvans fc Shmeld 


20131) . These authors describe models with a single 


mantle layer (one sublayer), only. This causes the model to be inaccurate. A major problem 
for such a model is, e.g., the interaction of CO, supposed to be in the outer layers, with the 
photodissociation products of H 2 O in the inner layers. The inclusion of several sublayers 
removes this and similar discrepancies, and can be considered as a novelty. Additionally, the 
sublayer approach allows to perform a limited analysis on the depth-dependent composition 
of the ice layer (section I3.5.3p . 

The rate of chemical reactions and species’ interchange between the layers is largely 
governed by molecular diffusion. The diffusion rate is dependent on the binding en e rgy E h 


of each molecule. In the case of surface species, the approach by 


Garrod fc Herbst 


( 120061 ) 
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was used, i.e., Eb^s = 0-5 Ed- 

For subsurface species, several differences have to be taken into account. First, unlike 
surface species, molecules in the bulk are surrounded by other molecules, and cannot 
directly desorb into the gas. This means that the term ‘desorption’ energy is fictional; in 
fact, it is the absorption energy Eb- By definition, Eb or Eb is the energy required to 
break all physical bonds for a molecule on ice, or within ice lattice. It is used to calculate 
the species’ vibrational frequency, z/q, Eq. ([ 8 ]). 

) assumed that Ejj or Eb are equal. However, the number of neighboring 
species in the lattice is higher than that for surface species, and Eb can be expected to be 
correspondingly higher than Eb- The precise Eb/Eb ratio probably should depend on the 
average number of species surrounding a molecule on surface and in bulk ice. This depends 
on small-scale surface roughness and ice porosity. Additionally, the molecular composition 
of ice is important, because strongly polar species (e.g., H 2 O, HCOOH, NH 3 ) are able to 
form much stronger intermolecular bonds than molecules with a low polarity (e.g., CO, 
CO2, N2). 


GarrodI (^2ni3a 


A constant each species’ Eb for the su rface layer and Er for each of the sublayers was 


adopted, in line with previous studies (e.g. 


Haseeawa et ah 


199211 . It was assumed that 




= S.OEn and E 


B,M'i 


= E 


B,MS 


= 3.3Eb for sublayers 1 , 2 , and 3, respectively. Eb 


for the inner sublayers has been assumed higher than that of sublayer 1 , because a major 
part of these sublayers consist of the polar H 2 O molecules, while sublayer 1 has a high 
proportion of CO after the freeze-out process has ended. In the present study, the difference 
between Eb,mi and Eb^m 2 i Eb^ms was taken to be mediocre 10%. Such a cautious approach 
was chosen to ensure that chemical reactions are possible in all ice layers, corresponding 


to photoprocessing experiments with astrophysic al ice analogs (e.g . 


For mantle binding energy E^^m the approach of 


Garrod &: Herbst 


Gerakines et ah 


3 


199611 . 


fl2006l) was retained, i.e.. 




















Eb,M = 0.5 Eb. The high Eb values imply a very slow diffusion, the molecules are basically 
frozen in place at temperatures relevant for this model fsection l2.5.3p . 


In addition to the above discussion, binding energy of mantle species must be higher 
than desorption energy of surface species, i.e. E^^m > Eb- This is to ensure that the 
movement of bulk ice species is slower than evaporation, a neces sary condition for the ice 


to be described as a solid. This condition was not recognized by 


GarrodI f)2013al) . In a case 


with Eb^M < Eb, ice molecules become mobile before evaporation, forming a liquid, which 
is not physically justifiable for interstellar or circumstellar conditions. 

The desorption energy for surface species can be dependent on the fraction of the 
surface consisting of H 2 molecules. Species, adjacent to a H 2 -covered surface may have much 
lower actual Eb than species attached directly to ice surface. In the present stud y, this 


was att ributed to the light H and H 2 species, only, with the approach proposed by 


Garrod 


fj2013al) . It was not attributed to heavy-atom containing species. In other words, heavy 
molecules do not ’step on to’ H 2 when accretion or diffusion occurs. Instead, the mobile 
hydrogen species make room and allow the heavy molecule to be always firmly bound to 
the ice surface. The vibration frequency of a heavy species, adjacent to H 2 , is orders of 
magnitude higher than the diffusion rate of H 2 . When the adjacent H 2 molecule hops away, 
no other H 2 molecules will be able to move into the now free adsorption site before the 
adjacent heavy molecule approaches the site via vibration, and binds itself strongly to the 
surface. 


2.4-2. The formation of the ice layer 

Figured] shows the ice layer, as it is described in the model. The formation of ice starts 
with the accretion of species onto grains. This forms the surface layer. If the thickness 
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Fig. 1.— Schematic drawing of the ice mantle on an interstellar grain as considered in 
the model. The allowed molecule transfers between the four ice layers are shown. In ice 
formation (‘compaction’) process, molecules are sequentially transferred from the gas to the 
surface S and, then, the sublayers Ml, M2, and M3, as long as net accretion happens. 
The molecules start to ‘pile up’ in the inner sublayers, up to a thickness of 60ML. The bi¬ 
directional diffusion of ice species is permitted at all times between adjacent ice layers. (A 
color version of this figure is available in the online journal.) 
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of the surface layer exceeds 2ML, the molecules are transferred to sublayer 1. They are 
then transferred to sublayer 2 and, after that, to sublayer 3, if their thickness does not 
exceed 60ML. This means that, initially, the molecules are transferred through sublayers 1 
and 2, and pile up in sublayer 3. After the abundance of molecules in M3 has reached the 
equivalent of 60ML, sublayer 2 is hlled up to the same thickness. All the sublayers reach 
an approximately similar thickness after the accretion period has ended (hnal ice thickness 
160-170ML, depending on the particular desorption model). 


The transition between the surface layer and sublayer 1 is of particular concern, 
because the chemical processes differ greatly in these two phases. Molecules in the surface 
layer are able to freely diffuse around, while they are mostly frozen in place in the case o f 


mant l e layers (section 12.Sp . Diffusion experiments with interstellar ice analogs fjPalumbo 


20061 : 


Rant et af 


2007aJ) have shown that molecules are highly mobile in porous ice layers. 


i.e., the surface area of ice is higher than that of a single smooth monolayer on top of the 
ice. This means that the nominal thickness of the surface layer is not limited to a single 
monolayer. 


2.4-3. Ice compaction: transition of species from surface to mantle 


Experiments (jPafumbp 


2006; 


Rant et ah 


2007b 


Si; 


Accolla et ah 


201 ll) have shown 


that interstellar ice analogs are porous, when obtained by depositing molecules onto cold 


surfaces. Ices obtained by c 


have a compact structure (jOba et ah 


eposition 


of free atoms and subsequent reactions on the surface. 


20091) . Additionally, it has been demonstrated by 


detailed models that slowly accumu 
accreted ice can be a highly porous I 


ating ices hav e a compact structure, while rapidly 


Garrod 


2013bl h 


Based on these evidences, it is possible to make out two ‘modes’ for interstellar ice 
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structure. The first is water-dominated ice, which forms by the slow accretion of H, O and 
other atoms, followed by reactions on the surface. This mode dominates during the long 
low-density period (uh < 2 x lO^crn”^ in core contraction. The second mode sets in as the 
contraction is accelerated and gas density exceeds lO^crn”^. The majority of species being 
accreted are CO and N 2 molecules that simply accumulate onto the grains with no chemical 
reactions synthesizing major ice (surface) components. It can be seen that the hrst mode 
corresponds to compact ice, and the second mode - to porous ice layer. The the formation 
of the subsurface ice has been described in the model with the intention to describe both 
modes, and the transition between them. 


It has been demonstrated that interstell ar ice analogs undergo compaction (reduction 


of p orosity), when 


10 ns 


Palumbo 


2006 


e xposed to UV pho tons (jPalumbo et ah 


2010 


Rant et ah 


2007b|) . or exothermic reactions (lAccolla et al.N201lh . 


cosmic rays 


(fast 


Additio nally, observational e vidence suggests that actual interstellar ices are compact in 


nature flPalumbo et ah 


20101). Because of this, the nominal thickness of the surface layer 


(representing porous ice) was never allowed to become thicker than lOML. In the model, 
the transition of species from the surface layer to the compact sublayer 1 was described as 
a hrst-order process. The rate coefficient was found empirically: 


fccomp = 5 X 10 \ 


( 2 ) 


where Bs is the thickness of the surface layer in MLs. The rate coefficients for other 
transitions associated with ice formation in the sequence S' —> Ml —)■ M2 —>■ M3 were 
calculated with Eq. (|2]), too. 


The compaction process is initiat ed only, when surface ice t 


This is consistent with the Ending by 


lickness exceeds 2ML. 
Va.svunin fc HerbstI (j2013al) that more than one 


reactive surface monolayer is necessary to properly describe surface chemistry. Based on 
pre-calculated ice maximum thickness, B < 180ML, it was assumed that each sublayer has 
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a maximum thickness of 60ML. Because the inner layers are hlled hrst, sublayer 1 does not 
get full in the models described here. The result of the mantle formation process can be 
seen in hgures [2] and dH 


Hydrogen diffusion 


Most of ice species diffuse between the ice layers S, Ml, M2, and M3 via thermal 
hopping 1 section r2.5.3p . The diffusion of hydrogen species H 2 and H was treated separately 


- with the use of data derived 
line with the earlier studies by 


rom experiments and dec 


Kalvans fc Shmeldl ( 2010 


icatec 


201311 ■ 


models. This approach is in 


The hydrogen diffusion coefficient in ice D (cm^s can be expressed as 

D = (3) 


where E^ct (K) is the activation energy for diffusion (IStrauss et ah 


and for Ho were calculated from c 


(l2005l) . as specified in 


ata supplied by 


Strauss et al 


199411. T he pa rameters Do 


(1199411 and 


A wad et al. 


KalvansI (l2013[l . The adopted parameters are Dq h = 10 
Do,h 2 = 4-76 X lO'^cm^s"^ i7act,H = 185, and E„,r±.Hn = 90K. 


The temperature-dependent va’ 
coefficient, according to Eq.(19) of 


act,/f2 

ue of D is then used to calculate the diffusion rate 


Kalvans fc Shmeldl (1201311 . The time-dependent 


thickness of the sublayers and the surface layer was taken into account for this 
calculation. The permitted transitions between ice layers for hydrogen diffusion are 
S <—)■ Ml —)■ M2 t —> M3. 

Because the diffusion of H 2 is relatively fast, the maximum proportion of H 2 in ice has 


to be regulated . It can be reasona 


with in the ice (IStrauss et al 


5% (lOberg et al. 


Dly assumed that most of the hydrogen resides in cavities 


1994 ( 1 . The estimated porosity of interstellar ices is around 


2011a|l . Following this, it was assumed that the diffusion of H 2 towards a 
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sublayer stops, if the proportion of H 2 in that sublayer exceeds 5%. 


2.5. Ice chemistry 


2.5.1. Surface binary reactions 


Reactions on grai n and ice surfaces have been considered with the approach of 


Hasegawa et al. 


(IGarrod fc Panlv 


992 


The reaction rate was adjusted by re action-diffusion co mpetition 


20111 ) and modified rate-equations method (jCaselli et alJll998h . The rate 


coefficient for binary surface reactions (cm^s ^), is 


kij Gact (y ) (Rdiffji T -Rdiffj); 


(4) 


where is 17act the probability for the reaction to occur, related to the activation energy 
barrier E^. The rate, (s“^) at which molecule i scans the whole grain via thermal hopping 


Rdiffji ^hop,i/-^S' 


(5) 


where 


fchop,i = z/o,iexp(-Efe,i/T). (6) 

is the rate (s“^) with which molecule i migrates from one surface site to an adjacent one. 


Following iGarrod fc Paulyl (120 llh . in the case of reaction-diffusion competition, the 
reaction probability is 


h'act(G) = 


vniij) 


VK{ij) + k hop,i T ^hopj 


(7) 


where u = z/o,i + and n{ii) = exp(—Reactive desorption (IGarrod et ahl 120071 ) is 


applied by assuming that one per cent of all reactions products transit to the gas phase. The 
modihed rate-equation method has been app lied to all surface reactions, as implemented in 


the ‘ALCHEMIC’ code ( Semenov et al. 


2010h . 
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The synthesis of CO9 on very cold 


of a CO molecnle flGarrod fc Paulv 


<8K) grains via 0+H—>0H reaction on top 


201 ll) has not been included in the model. It is a 


specihc process, which has not been attributed to other species and thus may introduce 
bias in calculation results. Additionally, as recognized by the same authors, molecules on 
surface can be adjacent to several neighbors. The O+H+CO mechanism requires that an 
O atom and, half-way into the reaction, an OH radical has attached itself exclusively to 
a single CO molecule. Other common interactions for the newly formed and still excited 
radical could include reaction with H2, which is almost omnipresent on cold grains, strong 
hydrogen-bonding or reaction with the abundant H2O, or even a reaction with CO2. The 
O-I-H-I-CO reaction has been veri hed as effective by the microscopic Monte-Carlo modeling 


technique (jChang fc Herbst 


20121) . but, again, their model by design implies that a surface 


molecule is bound to a single species below it. If OH and CO fail to react immediately at 
8K, this may well mean that OH becomes hxed in a conhguration in which it is unable to 
reach CO before it is tran sformed into H2O by H atoms. H has a non-negligible abundance 


even in dense cloud cores (jColdsmith fc Li 


200; 


). Based on these considerations, I suggest 


that it is acceptable not to use the tertiary reaction for CO2 synthesis in the current model. 


2.5.2. Chemical reactions in sublayers 


Due to the relatively high Eh^M values (> l.bEu, section [2.4.ip . molecules in bulk ice 
are practically immobile in the temperature range considered. Even light radicals with 
Ed = 800K require > 10^^ years to make a single hop at 16K, the maximum temperature 
permitted in the model. Because of this, one cannot assume that subsurface molecules 


and rac 


Garrod 


i cals ar e mobile within their respective sublayer, which is the case of the model by 


fj2013a|) . A new approach is proposed here, taking into account that molecules are 


locked in their absorption sites for long timescales. 
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Species in the mantle were treated as particles adsorbed to the surface of ice lattice 
cells. Because of the stronger binding of molecules in bulk ice, the characteristic vibration 
frequency z/qj fo r molecules in su b layer / is different from that of surface species. Following 


the formalism of 


Hasegawa et ah 


fll992h : 


z/Qj = 


( 8 ) 


where is the density of adsorption sites in a lattice cell, m is the mass of the adsorbed 
particle, and / is either Ml for sublayer 1, M2 for sublayer 2, or M3 for sublayer 3. 
Molecules in ice lattice cells vibrate with frequency uqj. A chemical reaction may occur, 
if they approach another molecule in the cell and overcome a certain energy barrier. 


E, 


prox* 


This ‘proximity barrier’ was assumed to arise from the adsorption force exerted by 


other neighboring species to the molecule that is approaching its reaction partner. Such a 
parameter has never appeared in the astronomical literature before. It cannot be chosen 
arbitrary, because it affects the abundance of chemical radicals in ice. It was assumed 
that F^prox = which results in radical (mainly OH and NH 2 ) abundances < 0.2% of 

total ice. The factor 0.1, multiplied by species’ sorptio n energy in ice, was us ed also for 


calculating lateral bond strength between molecules by 


Chang &: Herbst 


(1201211 


The scanning rate of the full surface of a single lattice cell for a molecule i is 


Ec,f,i ^hop,f,i/E^c ^0,/,i®^P( F/pj-ox/F')/A^,., (9) 

where Nc = 10 is the assumed average number of neighbors available for reactions in each 
lattice cell. The rate coefficient for a binary reaction (cm^s“^), in ice mantle lattice cell is 

^f,ij Fact,/(0) T (10) 

The reaction probability U^ctj is calculated similarly as in the case of surface reactions, 
Eq. ([7]). The modihed rate equations method was not applied in the case of mantle 
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chemistry because there is no accretion and no desorption in the bulk ice, new species arrive 
mostly via the photoprocess, and, in terms of this model, the reactive volume consists of 
a large number of interconnected lattice cells. Such a regime has not been tested by the 
Monte-Carlo approach. 


2.5.3. Diffusion of species in ice mantles 


Molecule thermal diffusion between adjacent sublayers, negligible for tempe ratures in 


the pr esent study, can be of high importance for ices in protostellar envelopes flGarrod 


2013al) . and thus was included in the model. The permitted directions for inter-sublayer 
diffusion are S i—> Ml <—y M2 <—> M3. Diffusion rate coefficients are calculated 
according to 


kd^ff — 


uojoexp{-Ebjo/T) 

yBfo 


( 11 ) 


where BfQ is the actual thickness (in MLs) of the source layer (S', Ml, M2, or M3), and y 
is the probability for a molecule to move one ML towards the target layer. In other words, 
y is the number of steps required for a molecule to move closer by one ML towards its 
target layer. Assuming a simple, cubic symmetry we get y = 6. The exact value of y is 
insignihcant, because bulk diffusion is very slow. The above means that yB^Q is the total 
number of steps for a molecule for the diffusion into an adjacent ice layer. 

Eb for molecules in the surface layer is much lower than that of mantle species. Because 
of this, surface species would diffuse much faster from surface to sublayer 1 than in the 
opposite direction. This is unrealistic, because the diffusion of molecules between ice layers 
involves ‘making room’ in the target layer (molecule swapping. 


Favolle et ah 


2011al) . Taking 


this into consideration and to reflect the diffusion process more precisely, I assumed that in 
the case of surface molecules, Ebjo = 7^6,mi for equation (fTTj) . This ensures an similar, slow 
diffusion rate between for the directions S —> Ml and Ml —> S. 
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2.5.4- Dissociation of ice species 


Dissociation by interstellar and cosmic-ray-induced photons has been included for gas, 
surface, and sublayers. For species in the icy mantles below the surface, the attenuation of 
radiation by overlying ice layers has to be considered. Each sublayer has a different (and 
time-dependent) number of overlying monolayers. For surface molecules, the attenuation 
factor was assumed unity. For subsurface species, the attenuation factor was calculated for 
the middle monolayer in the respective sublayer for each time-step. It was then attributed 
to all species in the particular sublayer. Taking into account that each monolayer has an 


absorption probability of Pabs = 0.007 flAndersson fc van Dishoeck 
factor is 

4ce = (1 - 


20081) . the attenuation 
( 12 ) 


where Pa is the number of monolayers above the sublayer in consideration, and P/ is the 
number of monolayers in this sublayer. 

Dissociation products remain in their parent sublayer. It was assumed that the 
recombination of the dissociated fragments is not a particularly preferred pathway, and is 


j ust a s possible as any other reaction. This is in line with the Endings by 


Andersson et ah 


fl2006l) that some of the fragments of photodissociated molecules tend to travel a few 
molecules’ worth away from their site of origin. Thus, the recombination possibility for the 
fragments is similar to the possibility for reactions with other species in that sublayer. 


Results 


As a benchmark reference sheet for the evaluation of calculation results the observational 
results by 


Whittet et al. 


(120071 ) are used. They provide the relative proportions of the 


major ice components - H 2 O, CO, and CO 2 - in samples of interstellar molecular gas 
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towards background stars in Taurus complex of dark clouds. Because the observational 
studies provide a visual extinction value for each held star, the results are presented in an 
y4y-dependent manner, too. The objects are listed in Table [TJ Each of them has an unique 
Ay, and I use their Ay values as an identiher for subsequent results tables. For additional 
result evaluation, the ice component abundances at threshold Ay - specihed in Table [U - 
are used. 


The use of data from 


Whiffet et af 


(120071) allows an evaluation of calculated ice 


composition that is more 

detailed than that c 

)f other recem 

; dark core models. For example. 

Garrod & Paulv 

(2011 

); 

Vasvunin & Herbst 

(2013a 

); and 

Chang & Herbst 

(2014 

) use a 


comparison at arbitrary 
observational results by 


points in time of the simu 


Booeert et ah 


(1201111 and 


ation run wit 


Oberg et ah 


1 Elias 16 or the averaged 


(l2011bl) . Contrary, the 


observations of ices towards eight background stars have been used for comparison in the 
present study, and the Ay of each star was used to tie the observational data to a particular 
point in time in the model. 


The methods and results for the different selective desorption mechanisms were 
investigated on a case-by-case basis. The Standard model, described in section |2l was used 
as a matrix, where a single parameter was changed and its effects evaluated. For an easier 
perception of the meaning of model results, results tables present the calculated-to-observed 
abundance ratio. The C0:H20 and C02:H20 ratios are considered (see note 2 below 
Table [T]) . 


3.1. The Standard model 

Figure [2] shows the calculated the evolution of ice structure. A slow ice accumulation 
starts from the earliest stages of the model because the low photodesorption yield is 
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Table 1: The observational data - Ay and abnndance ratios - used for compari son with calcu¬ 
l ation results, and the corresponding results of the Standarc 


(120071 ) ■ Table 1, unless otherwise noted. See 


Whittet et ah 


model. Data from 


Whittet et ah 


(2007) for more details. The cal¬ 


culated abundance of H 2 O and CO 2 in monolayers for their respective Ah is given in the 


lower part of the table. 






Observations 

Standard model 

Source ID 

Association 

Ay 

Time^ , kyr 

CO 0 / 
H 2 O v® 

CO 2 0 / 
H 20 A 0 

calc 2 
obs 

C02^ 

043728.2-h261024 

Tamura 2 

6.3 ± 1.5 

932 


25.0 


3.8 

042324.6-h250009 

Elias 3 

10 ± 0.5 

958 

20.2 

19.1 

1.7 

3.1 

043325.9-h261534 

Elias 13 

11.7 ± 0.5 

964 

12.0 

16.7 

3.2 

3.1 

043926.9-h255259 

Elias 15 

15.3 ± 0.5 

973 

27.3 

16.7 

1.7 

2.6 

042630.7-h243637 


17.8 ± 1.5 

976 

45.1 

18.3 

1.1 

2.2 

043213.2-h242910 


20.9 ± 1.5 

980 

33.3 

16.0 

1.6 

2.4 

044057.5-h255413 

Tamura 8 

21.5 ± 0.5 

981 

24.3 


2.2 


043938.9-h261125 

Elias 16 

24.1 ± 0.5 

983 

25.3 

21.0 

2.2 

1.7 






H 2 O, ML 

CO 2 , ML 

CO, ML 

H 2 O threshold^ 


3.2 ± 0.1 

838 


6.7 



CO 2 threshold 


4.3 ± 1.0 

894 



12.0 


CO threshold^ 


6.8 ± 1.6 

938 




3.5 


^Time required to achieve the particular Ay in the model 
"E qual to nlnhs 


Whittet et al 


(l2nnih 
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Fig. 2.— Ice accumulation (total thickness in monolayers) on grains during late cloud evo¬ 
lution for the Standard model. S is the surface layer, Ml, M2, and M3 are sublayers 1, 2, 
and 3, respectively. 














- 33 - 



-HCOOH 

= = =H2C30 
-CH3CHO 

HCOOCH3 

-CH2CO 

-CH30CH3 

.CH30 

.HCO 

time, Myr 



-HCOOH 

= = =H2C30 

-CH3CHO 

HCOOCH3 

-CH2CO 

-CH30CH3 

.CH30 

.HCO 


time, Myr 


Fig. 4.— Gas-phase (top panel) and ice (bottom panel) abundances, relative to hydrogen, 
of the Standard model for organic species observed in quiescent dark clouds. The ‘tremble’ 
of some abundance curves at early times is an artifact, related to the step-like transition of 
surface species into sublayer 1. 
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insufficient to keep the grains clean at Ay ^ 2 . 

Table m shows that both, CO and CO 2 , are signihcantly overproduced. Such a result 
was obtained, partially because the calculation results are compared to observational results 
tied to a particular Ay, and not to single averaged values, which is the case in many 
previous papers. This implies the assumption that the observational sample consists of 
low-mass prestellar cores. They are in different stages of evolution, or, alternatively, the 
held stars sample molecular gas in different distances from the center of the core. 

The thickness at which photodesorption achieves its maximum efficiency (i?ph) is 
approximately 3MLs in this model. Test calculations with different Yp^ reveal that a grain 
can be covered in ice with thickness lower or similar to Sph for very long timescales at Ay 
lower than the water threshold value. 


Water ice abundance of seven monolayers at A^y, is signihcantly more t 


ice layer adopted as a limit for water ice detection by 


ran the single 


Garrod fc Paulvl f 201ll) . The total ice 


thickness at 3.2mag is 27ML, with HCOOH, CO 2 , N 2 , and NH 3 being ot her major spec ies. 


’his result of the Standard model can be consistent with observations. 


Whittet et ah 


( 1200 ll) conclude that water At^ is coincident with an increase of grain size, most probably 


by ice mantle accumulation. Approximately 20 to 30 monolayers of ice are require c 


moderate increase of the grain radius by ten per cent. In addition. 


Whittet et ah 


for a 


(1200 ih 


allows the possibility that the ice is concentrated in a clump along the line of sight. This 
would imply that the actual Ath is halved, although such an extreme case is unlikely. Both 
of these considerations permit the existence of water ice at lower Ay or thicker ice at 
Ath. Water ice abundance, relati ve to hydrogen, reaches 10"*^, which is in agreemen t with 


observations of molecular clouds (IBoogert fc Ehrenfrennd 


2004; 


Whittet et ah 


2003 ). 


Figure |3] shows the evolution of the abundances o 


with observational data flGibb et al 


2004; 


Oberg et al 


major ice species. A comparison 


2011a . Elias 16 with Ay 24.1) 
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reveals that ammonia is overproduced (NH 3 :H 20 = 21 %), while methanol is underproduced 
(CH 30 H;H 20 = 0 . 1 %). NH 3 is synthesized by a surface reaction pathway, similar to that of 
water, and both compounds have very similar E^. Thus, ammonia is almost inseparable 
of watery ices in the model, and always has a relatively high abundance. Methanol has 
a low abundance largely because its synthesis is hampered by barriers in formaldehyde 
hydrogenation reactions. These effects are unlikely to significantly affect the proportions of 
H 2 O, CO, and CO 2 , although they are important for the chemistry of nitrogen and complex 
organic molecules (COM). 

Early ice accumulation at relatively low Ay results in a phenomenon that can be 
termed as ‘photon-dominated ice’ (PDI). This means that the (relatively thin) ice layer, 
formed before a high extinction of the ISRF is reached, experiences intensive irradiation 
by interstellar photons. PDIs are promoted by low yields for Ay-dependent desorption 
mechanisms - interstellar photons and reactive desorption. The Standard model (along with 
the low-yield reactive desorption Model B) is among the most pronounced representations 
of PDIs described in this paper. 


In the model presented here, photoprocessing does not affect significantly the 
abundances of major species because radicals in ice regenerate these species. However, 
the production of organic species is greatly enhanced in bulk PDIs. Ice photoprocessing, 
coupled to non-thermal desorption mechanisms can be a major source of gas-phase COMs 
in quiescent clouds. During the long-lived diffuse phase (uh ~ 3000cm“^) it produces 
abundances on the order of 10 ~^ ^... 10 —relative to hydrogen for species that recently hav e 


been observed in quiescent gas fjOberg et ah 


2010 


Bacmann et ah 


2012 i 


2012 


Cernicharo et ah 


, see also Table | 8 ]). Figure 0] shows the calculated abundances for several of these 


species. 


The gaseous COM-rich evolutionary phase has a length of almost IMyr in the Standard 
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model. It starts with the formation of the hrst sublayer at ISkyr and ends at densities higher 
than 10^cm“^. After this period, the rate of species’ transition into sublayers is faster than 
the synthesis of COMs, which takes several steps to be completed. COMs become buried 
in ice, and are not available for desorption in the near-surface layers. This mechanism 
inhibits desorption for any complex species during the freeze-out stage, regardless if they 
are produced in the surface or the mantle. 


The Standard model rather poorly represents ice accumulation, and the COM-rich 
period is signihcantly shorter for models with effective desorption mechanisms. Because 
of this, a detailed discussion on the chemistry of COMs has been reserved for the 
‘complete’ model fsection l3.5.ip . which describes overall ice composition more accurately. 
Earlier papers th a t inve s tigate the chemistry and reaction ne t works of COMs include 


Garrod fc Her 


Taquet et al 


3St 


fl2nn6h: 


(12Q121); 


Garrod et al 


Vasvunin fc Herbst 


1200811: 


Belloche et 


fcoiObfl : 


al 


1120091); 


Laas et al. 


(|201l|); 


GarrodI fl2013al) . The gas-phase organic- 


rich period has little effect on most carbon-chain species because they are not produced in 
the sublayers. 


3.2. Photodesorption 


Desorption by ISRF and secondary photons ha. s been often treatec 


mann e r in recent papers flVasvnnin fc Herbst 


2013a 


bl; 


Chang &: Herbst 


in a 

non-selective 

2014; 

Gerner et ah 


20141) . iGarrod fc Paulvl fj201ll l use experimental photodesorption data, obtained for pure 


species (water, carbon oxides, nitrogen), in their models, which help to explain Ath for 
different ice constituents. These authors adopt similar desorption yields for ISRF and 
secondary photons. 


Experiments have shown that codesorption from icy mixtures may result in similar 
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Table 2: Photodesorption yields as derived from experimental data. 



fhsrf 

^rph 

H 20 

2.7E-03 

1.8E-03^ 


5.7E-03 

3.0E-03 

C 02 

3.5E-03 

2.3E-03" 

N 2 '’ 

5.5E-03 

3.0E-03 

02 " 

3.3E-03 

2.6E-03 

Other 

3.0E-03 

2.0E-03 


" Oberff et al. 

(2nn9al 

^ Bertin et al. 

(201,1) 

' Obere: et al. 

(2nn9bi 

Favolle et al. 

(201,11 
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Time, Myr 

Fig. 5.— Calculated abundances of major ice species, relative to hydrogen, photodesorption 
simulation with empiric, i?£)-dependent yield. Compare with £gure[3l 
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Table 3: Comparison of observations and calculations with modified photodesorption yields 
for ice species: abundances of ice species CO and CO 2 , relative to water, and abundance in 
monolayers for H 2 O and CO 2 at threshold Ay- Photodesorption model A with uniform yield 
for all species. Models that use experimentally detected photodesorption yields or empiric, 
A’^-dependent yields give very similar results (see text). 


A 


Ay 

CO£^ 

obs 


6.3 


3.9 

10 

1.9 

2.8 

11.7 

3.6 

2.8 

15.3 

1.8 

2.3 

17.8 

1.2 

2 

20.9 

1.7 

2 

21.5 

2.3 


24.1 

2.3 

1.5 

Ath 

H20:C0:C02(ML) 

3.9:3.9:10.2 
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yield for different species flOberg et al.ll2009bl: iBertin et alJl2013ll. and that 


efficiency is wavelength-dependent flFavolle et al. 


2011 


2013 


Berlin et al. 


photo desorpt ion 


2013 ). In order 


to investigate these factors, three simnlations have been rnn, each with a different approach 
on photodesorption. 


The first simulatio n made use of experi ment al data, wher e 


desorption the dat a by 


authors, as well as 


Favolle e 


Oberg et al 


; al 


fl2009al 


l|2qi3h and 


Berlin et ah 


possi ble. For ISRF-induced 


(120131) was used. These 


bl), also provide data for desorption by secondary 


photons. Molecules with available yields are H 2 O, CO, CO 2 , N 2 , and O 2 . I could not 
hnd recent information on Risi-f for H 2 O and CO 2 . From the experimental data, it can be 
estimated that Fisrf ~ l.Sycrph- For species with no data available, the yields have been 
estimated. It was assumed that they are slightly higher than water photodesorption yields, 
and retain the 1.5 ratio. These data are summarized in Table O All experimental results 
are relevant for a minimum temperature of 15K and, and, for consistency, all the yields 
used in this model are those for 15K. The yields also assume mantles of inhnite thickness. 
Some thickness-dependent ef fects have been alr eady included into the model (section l2.3p . 


Bertin et ah! (120131) take into account the codesorption effect 


Only the yields provided by 
(for a C 0 :N 2 mixture). 

In the second simulation uniform photodesorption yields for all species were assumed 
(Tabled last row), i.e., codesorption is prevalent. 

For the third simulation, a weak, empiric yield dependence on Eo was used: 

blsrf = Qoexp(-F;i)/(Qi-MOOT)), (13) 

where Qq = 0.007 and Qi = 4200. Eq. (fT3ll approximately matches the experimental 


at a and, by 


(lOberg et al 


extrap olation, fits the temperature dependence for water photodesorption 


2009b|) in the relevant temperature interval (6-16K). It was applied for all 


surface species. As in the hrst two simulations, Idsrf = l.STcrph- 
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All of these photodesorption models produce very similar abundance patterns, and an 
example is shown in hgure O The abundance of carbon oxides in ice, relative to water 
at different Ay values (Tabled]), is very similar for all three simulations, as well as the 
abundance of these species at Ah- These numbers do not deviate by more than 5% for 
the three simulations. Because of this, the results of the uniform-yield simulation (Model 
A, see below) have been presented in Table [31 only. The effective desorption yields in real 
interstellar conditions should be much more selective than experiments indicate in order 
to introduce signihcant changes for the H20:C02 ratio. Based on current evidence, it can 
be concluded that, although photodesorption largely regulates molecule accretion on grains 
under diffuse-cloud conditions, it can be regarded as non-selective because it basically does 
not affect the relative proportions of ice ingredients. Similarly to the Standard model, 
models with modihed photodesorption yields show a period of a relatively high COM 
gaseous abundances. Because of a higher Tph, this phase begins much later, at around 
0.6Myr. 

Based on these results, I conclude that it is unnecessary to use a complex, selective 
approach for molecule photodesorption from interstellar grains. Thus, the simplest approach 
with uniform desorption yields (0.002 for secondary and 0.003 for ISRF photons) was chosen 
as the most appropriate photodesorption model. This is further referenced to as Model A. 


3.3. Reactive desorption 


Reactive desorption was included in the model following the approach by 


Garrod et al 


(120071 ). based on the Rice-Ramsperger-Kessel (RRK) theory. For clarity, a full description 


is included here. The fraction of reactions resulting in product desorption is 


a-pRRK 
1 -I- oPrrk ’ 


( 14 ) 
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Table 4: Comparison of calculation results. Reactive desorption models B (q;=0.01), C 
(q;=0.03), and D (q;=0.06). 




B 


C 


D 

dly 

obs 


obs 


obs 


6.3 


3.8 


4.0 


4.2 

10 

1.7 

3.1 

1.8 

3.1 

2.0 

3.2 

11.7 

3.2 

3.1 

3.5 

3.2 

3.7 

3.2 

15.3 

1.7 

2.6 

1.8 

2.6 

1.9 

2.6 

17.8 

1.1 

2.2 

1.1 

2.2 

1.2 

2.2 

20.9 

1.6 

2.4 

1.7 

2.3 

1.8 

2.4 

21.5 

2.2 


2.3 


2.5 


24.1 

2.2 

1.7 

2.3 

1.7 

2.5 

1.7 


HaOiCOiCOs (ML) 

HsOiCOiCOa (ML) 

HaOiCOiCOa (ML) 

Ah 

6.9 : 

3.5 : 12.1 

5.6 : 

3.6 : 11.5 

4.5 : 

3.7 : 10.9 
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where a = z/q/i^s is “the ratio of the surfa cemolecule bond fre quency to the frequency at 


which energy is lost to the grain surface” flGarrod et ah 


2007 ). The probability Prrk for 


an energy, larger than Eo-, to be present in the molecule-surface bond depends on the total 
energy -Ereac released in an exothermic reaction: 

Ed 


-pRRK — I 1 — 


Er. 


S-1 


( 16 ) 


where s = 3N — 5 and N the number of atoms in the most complex product molecule (s = 2 
if At = 2). 


Ej-eac is equal to negative reaction standard enthalpy, AH^. In the model, the latter 
was calculated for each reaction, using the gas-phase standard enthalpies of formation AHj 
of reactants and products. Standard enthalpies were chosen, because more species have 
thermochemical data available for 273K than for 0 or lOK, while the changes are relatively 
insignihcant. Similarly, enthalpy, not Gibbs energy, was used because entropy is of little 
importance at the low temperatures considered in the model. 


In order to facilitate the implementation of the reaction-specific reactive desorption 
into the program, only two cases were considered - either that all products desorb into the 
gas or that all products remain on the surface. This is the original approach used in the 
ALGHEMIG code. Consequently, the sum of product desorption energies were used instead 
of Ed in Eq. ffTSjl . N used in the same equation was assumed to be the number of atoms in 
the most complex species among the products. Although not entirely physically justified, 
these approximations generally retain the unique selectivity of the reactive desorption 
process. Because of the dependence on two variables - E^^ac and Ed - reactive desorption 
is different from other mechanisms, which have Ed as their only variable, unique for each 
species. 


The use o 
assumption by 


actual Erpar and E n values for each reaction in Eq. flT^ shows that the 


Garrod et ah 


(120071 ) that Prrk ~ 1, and thus /rd ~ is seldom applicable 
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for surface reactions. If a is taken between 0.01 and 0.1, less than a quarter of all reactions 
have /rd > 0.5q!. Usually, these are oxidation-reduction reactions with high exothermicity 
between light species. Atom exchange reactions involving heavy, complex molecules tend 
to have insignificant desorption probabilities. The changes implemented in this model 
mean that reactive desorption is a truly selective deso rption mechanism, with a as its 


only poorly-known parameter (see 


Garrod et ah 


20071. for a more detailed discussion). 


Calculations were performed with a values in the range 0-0.2, and Table H] shows results for 
0 =0.01 (Model B), 0.03 (C), and 0.06 (D). 


It can be seen that in terms of carbon oxide abundances relative to water. Model B is 
the best choice, although it still has a significant overproduction of CO 2 . The low efficiency 
of reactive desorption in Model B means that the total rate of desorption is actually lower 
than in the Standard model. The higher a values for Models C and D result in an even 
higher C02:H20 ratio. All the formation reactions of these two molecules are highly 
exothermic and result in high desorption efficiency (/rd > 0.7a). Because water forms in two 
steps, while CO 2 forms in two parallel one-step reactions, water is desorbed more effectively 
than carbon dioxide. These results mean that the model presented here does not place new 
constraints on a, although lower values may favor a C 02 :Il 20 ratio that is slightly closer to 
observations, as shown in Table 01 CO is formed by gas-phase reactions and is practically 
unaffected by reactive desorption. 


Selective reactive desorption favors the formation of COMs in ice, with COM 
abundances of Model B being higher than those in the Standard model results. This is 
because the high Ed of COMs combined with often low Uj-eac result in typically very low 
desorption efficiency for most reactions involving organic species. However, at high a 
values the ice does not accumulate until late times, and the length of the COM-rich cloud 
evolutionary phase is greatly reduced, if a Ri 0.1. Such high a values are not supported by 
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current theoretical evidence flGarrod et al. 


20071). 


3.4. Desorption by H + H reaction on grains 


The rate coefficient for desorption of ice molecules following an exot 
reaction on grain surface (indirect reactive desorption, H 2 FD for short, 
is proportional to H 2 formation rate per grain per adsorption site: 

[H] 


lermic H 9 formation 


Willacy et al 


199411 


^H2fd — -Racc.H X /H2fd X 


^UgNs 


-1 


(16) 


where i?acc,H, is the rate of accretion for H atoms onto grains and /H 2 fd is an efficiency 
parameter - the number of desorbed molecules per formed H 2 molecule. The last term 
in Eq. ([TB]) is the quadratic abundance of gaseous H atoms per grain per adsorption site. 
Similarly to photodesorption and evaporation, H 2 FD can affect sublayer species, if the 
number of overlying MLs does not exceed two. According to Eq. f[T6|) . indirect reactive 
desorption considers any surface exothermic reaction involving H, with the formation of H 2 
being the dominant energy source for this mechanism. 

No experimental or quantum-chemical data on the possible value of /H 2 fd were found 
in the literatur e. The only es t imate s of fn^fd were found to be provided by astrochemical 


modeling, e.g.. 


Roberts et al. 


fj2007(l . These values are not u sed here. The physic al model 


Roberts et al. 


( 20071) and the 


and methods for H 2 FD rate calculation differ signihcantly in 
present study. 

For estimating for /H 2 fd, a best-£t method was used. This means that simulations 
with a range of fn^fH values w e re per formed until the closest possible agreement with the 


observations by 


Whittet et al. 


(120071) was achieved. 


For the se’ 


Roberts et al. 


ectiy ity of H 2 FD, the approach developed bv iDulev &: WilliamsI (119931) and 


20071) has been used. This means that all species are desorbed with an equal 
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efficiency /nafd, up to a certain adsorption energy threshold, i?th- This is probably the 
simplest possible form for a selective desorption mechanism. H 2 FD is poorly constrained, 
and I explore cases with three valnes, as well as empiric eqnations for /nafd- The latter 
have been developed taking into acconnt the resnlts obtained by calcnlations with the i?th 
approach. 


The three threshold cases are with En <1210 K, Ed <2600K, and Ed <6000K. The 


hrst value is the one used by 


Roberts et al. 


(120071) . For water and carbon oxides this 


means that only CO desorption is enabled in the hrst case, CO and CO 2 may desorb 
in the second case, and all three components desorb in the third case. The obtained 
calculated-to-observed ice relative abundance ratios CO and CO 2 are shown in Table [5l 
Additionally, two models with a more complex approach on i^^i-dependent desorption 
efficiency have been investigated. 


3.4-1- 121 OK desorption energy threshold 

For a threshold energy of 1210K, calculations with values in the range 10“® — 10“"^ 
were performed. The most likely value was found to be near 4 x 10“®. Results of Model E 
with /Hafd = 4 X 10“® are shown in Table [5l The model produces a rather good agreement 
with observations for the C0:H20 ratio. Although carbon dioxide is synthesized via CO, 
the C02:H20 ratio remains almost unchanged from in the Standard model. /Hjfd is roughly 
inversely proportional to the abundance of gas-phase atomic H. 


3.4-2. 2600K desorption energy threshold 


The ap plication of higher va lues of Eth is an extrapolation of the considerations 


discussed by 


Roberts et al. 


(I2nn7h and 


Dulev fc WilliamsI (1199311 . This is partially justihed 
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Table 5: Comparison of calculation results. Indirect desorption by H+H —)■ H 2 reaction for 
molecules with adsorption energy threshold approach: < 1210K and = 4 x 10“® 

(Model E), Ed < 2600K and fn^td = 2 x 10“® (Model F), and Ed < 6000K and fn^fd = 
2 X 10-® (Model G)._ 




E 


F 


G 

dly 

obs 


obs 


obs 


6.3 


4.6 


0.9 


1.5 

10 

0.5 

4.2 

0.9 

0.9 

2.9 

1 

11.7 

1.1 

4.3 

1.8 

0.9 

5.1 

1.1 

15.3 

0.8 

3.5 

1.1 

0.9 

2.3 

1 

17.8 

0.5 

2.9 

0.7 

0.8 

1.4 

0.9 

20.9 

0.9 

3.1 

1.2 

0.9 

2 

1 

21.5 

1.3 


1.7 


2.8 


24.1 

1.4 

2.2 

1.8 

0.7 

2.7 

0.7 


HaOiCOiCOs (ML) 

HsOiCOiCOa (ML) 

HaOiCOiCOa (ML) 

Ah 

6.4 ; 

1.1 ; 11.3 

7.1 : 

1.9 : 2.8 

0.5 : 

1.5 : 0.7 
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by the fact that lower values of fn^fd are required by the approach used here. A threshold 
Ef) < 2600K means that H 2 FD affects both, CO and CO 2 , and that they have equal 
desorption yields. With such an approach it is possible to achieve a much closer match for 
observation-modeling results of H 20 :C 02 :C 02 ice abundance ratio than for all previous 
methods. Namely, the calculated-to-observed relative abundance ratio varies only within 
a factor of two for CO and CO 2 . The C0:H20 ratio has a higher spread. Slightly lower 
efficiencies than those in the Eth = 1210K model are required to achieve the best agreement 
with observations. Table [5] shows the results of Model F with a 2600K threshold and the 
best-£t value for /nafd of 2 x 10“®. 

These results show that a H 2 FD mechanism that desorbs CO and CO 2 , and excludes 
water, can bring a solution for the Ay-dependence of the relative abundances for major ice 
constituents. However, a threshold of 2600K either underproduces CO 2 or overproduces CO 
ice. This behavior has the element of F^£)-dependence. In order to further investigate this, 
an empiric H 2 FD approach is presented in the next subsection. 

I would like to note that, similarly to all other models, CO ice abundance at Ath (in ML) 
is significantly lower to respective H 2 O and CO 2 values. Nevertheless, the Ay-dependent 
results middle column in Table [5] clearly show that CO is the most overproduced at Ath of 
the three major ice species. This suggests that the chosen Ath value might be inaccurate. 
We shall return to this issue in the Conclusions section. 


3.4-3. Other approaches on H 2 FD efficiency 

It is possible that H 2 FD affects molecules with an even higher Ed. Calculations were 
performed for a threshold energy of 6000K, which means that H 2 FD was attributed to 
water and ammonia. The calculations produced an excess of CO, CO 2 in an very thin ice 
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Table 6: Comparison of calculation results. Indirect reactive desorption efficiency calculated 
with the use of the RRK theory (Model H), and by assuming a thermal desorption approach 
(Model I). _ 




H 


I 

dlv 

obs 


obs 


6.3 


1.2 


1.2 

10 

0.5 

1.1 

0.8 

1.0 

11.7 

1.0 

1.2 

1.5 

1.1 

15.3 

0.7 

1.1 

0.9 

1.0 

17.8 

0.5 

1.0 

0.6 

0.9 

20.9 

0.8 

1.1 

1.0 

1.0 

21.5 

1.2 


1.5 


24.1 

1.3 

0.8 

1.6 

0.7 

^th 

HsOiCOiCOs (ML) 

7.0 : 1.0 : 4.6 

HaOiCOiCOa (ML) 

3.9 : 1.3 ; 2.6 
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(Model G in Table [5]). 

Selectivity achieved by assuming a desorption energy threshold is a rather crude 
approximation. Two cases were considered, where the desorption efficiency is estimated 
in a more targeted manner. First, the RRK theory, which has been used in the case of 
direct reactive desorption fsection IXT]) . was applied also for indirect desorption (Model H). 
This was based on an assumption that indirect desorption occurs, when the newly formed 
and highly excited H 2 molecule transits energy to (or ‘kicks’) a nearby molecule. Taking 
into account the calculation results with energy thresholds, the calculation parameters are 
chosen so that they produce /nafd close to 2 x 10“® for CO and CO 2 and are negligible for 
H 2 O. This means that /nafd is calculated by Eqs. ffTTj) and flT^ with a = 6 x 10“® and 
-E'reac = 6 X lO^K. Only molecules with Ed < Ereac can be considered. These parameters 
have been applied for H 2 FD, which was included the Standard model, described in Sect O 


Another approach is to treat H 2 FD as evaporation, i.e., to assume that the energy 
transit from H 2 results in an immediate temperature increase for the nearby molecule 
(Model I). The desorption efficiency is calculated with an Arrhenius-type equation 
/H 2 fd = Q 2 GW{~Ed/{Q 3 )■ The parameter Q 2 = 7 x 10“® can be expressed as 

Q 2 = l^ot{Q3)P}l2id, (17) 

where ^(Qa), s, is the time spent for the nearby molecule at temperature Q 3 = 1.7 x lO^K, 
and PH 2 fd is the probability for the energy transfer from H^. Regardless of this, the approach 
used here can be viewed as an empirical derivation, only. Quantum effects play a significant 
role on the level of a single molecule. 


Table | 6 ] shows that models H and I produce a reasonably good £t to the observational 


results by 


Whittet et ah 


( 2007 1. All the calculated C 0 :H 20 and C 02 :H 20 ice abundance 


ratios fall within a ±50% margin of the observed values. This is better than Model F with 
Eth =2600K, where it was found to be impossible to achieve a similar degree of agreement 
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for CO and CO 2 simultaneously (Table |5]). Additionally, Models H and I produce a much 
greater similarity for H 2 O and CO 2 abundances at their respective threshold Ay- 


S.4-4- Summary for H 2 FD 

From the results of Models E through I (Tables [5] and [6]) it can be seen that the best 
agreement with observations for the H 2 FD mechanisms is achieved by models with (a) 
effective desorption of CO and CO 2 , (b) no or very low desorption efficiency for H 2 O, and 
(c) CO desorbed slightly more effectively than CO 2 . The sequence in C0>C02 >>H20 for 
desorption efficiency implies a dependence on the adsorption energy Eq. A semi-empirical 
(Model H) or empirical (Model I) approach produces an almost perfect fit to the 
observational H 20 :C 0 :C 02 ratio. 

It has to be emphasized that the rate of H 2 FD is highly dependent on the gas phase 
abundance of atomic hydrogen, which, in turn, is determined by the the chosen physical 
parameters of the model (gas density and Ay). Other assumed properties, such as the 
number of surface MLs or desorption from subsurface layers, can be of importance, too. 
The obtained /H 2 fd values are valid for a model with n(H) 0.01?7,(H2) during the long 
quiescent cloud phase with uh ~ 3000cm“^. /nafd is roughly inversely proportional to the 
relative abundance of H. An atomic hydrogen abundance of 0.001 relative to H 2 results in 
that a model with Eth = 2600K with /rd = 4 x 10“® is the model that produces the best 
agreement with observations. The parameters a and Q 2 for Models H and I, respectively, 
also scale inversely to n(H)^. The major finding is that a C0~C02 >>H20 sequence in 
desorption efficiency is required for the reproduction of observational results and that H 2 FD 
is a likely candidate mechanism, which could exhibit such a selectivity. 

A general trend for indirect reactive desorption models is that there is relatively more 
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CO 2 than CO at low Ay and vice versa at high Ay valnes (see Tables |5] and |6]) . This is 
a tentative trend that shows an inverse dependence on Bd. Cnrrent experimental resnlts 
snggest that CO has the higher desorption yield for ISRF photons (Table [3]), while direct 
reactive desorption significantly hampers the accnmnlation of water. Thns, this apparent 
discrepancy cannot be removed by adjnsting the efficiency of other desorption mechanisms 
examined in this paper. It is possible that the trend arises from discrepancies in the 
treatment of H 2 FD or gas-phase chemistry. 

There is a lack of experimental evidence or detailed qnantnm calcnlations on exothermic 
reactions resulting in the desorption of nearby molecules. Currently, it seems that the 
feasible /nafd values used in the model do not contradict existing views on reactive 
desorption. For typical atomic H abundances in dark cores, the efficiency for H 2 FD is 
two to four orders of magnitude lower than efficiency of direct reactive desorption. The 
signihcance of this mechanism lies in the sheer number of hydrogen atoms sticking to and 
meeting on grain surfaces. Because of the extremely low efficiency, this process can be 
hard to verify experimentally. An experiment with a beam of atomic H aimed to an inert 
cryogenic ice layer, e.g., N 2 , would be advisable. 


3.5. The complete model 

Nine derivations of the Standard model have been presented so far, each focusing on 
the efficiency of a single desorption mechanism with different parameters. I continue with 
combining the mechanisms in order to create a model that fully and adequately treats 
molecule desorption during the ice formation epoch in molecular cloud cores. 

Evidently, H 2 FD has the most signihcant effect on the relative abundances of major 
ice constituents. The application of H 2 FD (section 13.4p produce calculation results within 
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Fig. 6.— Ice accumulation on grains during late cloud evolution for the complete Model J 
(cf. figure [2]). 
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±50% margin of the observed values for CO and CO 2 relative abundances in ice. Its 
efficiency parameters have to be carefully chosen. Desorption by ISRF photons is relatively 
well constrained and largely regulates the ice accumulation and thickness of ices during the 
early stages. Reactive desorption has the general effect of favoring small molecules that 
form in highly exothermic reactions involving free atoms. Its efficiency is governed by a, a 
parameter which is known probably within an order of magnitude. 

A general result for all models is that, in terms of this paper, a higher desorption 
efficiency produce thinner ices at early stages and thicker ices (i.e. more ice molecules) 
at the end of the integration time. This is because subsurface photoprocessing efficiently 
produces multi-atom species ~ including COMs - in ice. The abundances of these molecules 
are lower, if ices accumulate later, and are exposed to ISRF photons for a shorter time-scale. 
This means that more atoms are concentrated into simple two- or three-atomic species, 
which means a higher number of molecules. The Standard model produces a hnal ice 
thickness of 160ML, while the combined models produce approximately 170ML. 

In order to establish a complete, functional model, calculations with various 
combinations of model cases A through I have been performed. In the complete model, all 
the three different approaches on photodesorption produce abundances of major species 
that differ only by hundredths. These are experimentally detected yields, uniform yields 
for all species, or an empirical relation mimicking the experimental results. The preferred 
choice here is the simplest case with an uniform desorption yield for all species (Model A). 

For reactive desorption, an increase of a in the range of 0.01-0.06 results in higher 
proportion of CO and lower proportion of CO 2 for all Ay values (Table H]). It has the 
beneficial effect of reducing the abundance of ammonia ice, which is overproduced in all 
models. All the abundance changes are within a few per cent at most. A more pronounced 
effect on H 2 O, CO, and CO 2 abundances at Ath can be observed (Table 0]). 0.03 and 0.01 
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(Garrod et ah 

2007; 

Garrod & Paulv 

2011) 


these values are appropriate for the current model. 


For the indirect reactive desorption mechanism an obvious choice is the empiric 
approach of Models H and I. A model with F^th =2600K produces a slightly worse match 
with observations fsection l3.4.4p . Other H 2 FD threshold models fail, because a similar 
desorption efficiency for CO and CO 2 and limited or no desorption for H 2 O is a requirement 
for reproduction the observed carbon-oxide-to-water proportions. Additionally, Models 
H and I are probably physically more understandable than the single energy threshold 
approach for models E to G. 

Table [7] shows the results of the three most important combined models. The chosen 
parameters are (1) uniform photodesorption yields (Tlsrf = 3 x 10“^ and Yc^ph = 2 x 10“^); 
(2) for reactive desorption a = 3 x 10“^; and, (3) three variants for indirect reactive 
desorption, Ed threshold of 2600K (Model J), H 2 FD efficiency calculated by the RRK 
theory (Model K), and H 2 FD calculated empirically as for a thermal process (Model L). 

Out of all these simulations. Model K yields a ±50% agreement with observations for 
relative proportions of carbon oxides and water. Model L produces a very similar result. 
Both these models also produce similar H 2 O and CO 2 abundances that differ by no more 
than 15% at their respective Ath- This can be regarded as an additional positive verihcation 
of results. Figs. [6] and [7] show the main results of Model L - the evolution of ice sublayer 
thickness and the abundances of major ice species, respectively. 

All the calculations described above are for a low-mass molecular core, whose Ay 
increases slowly. In order to provide a simple test for conditions relevant for massive cores, 
results are presented for Model M with a contraction time of 0.2Myr, i.e., hve times shorter 
than other simulations. The evolution of other physical parameters from start to end was 
retained as described in section 12.21 Model K set of desorption mechanisms was employed in 
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this simulation. As shown in Table [71 Model M produces up to two times lower abundances 
of carbon oxides relative to water, when compared to other versions of the complete model. 
If we were to obtain a better agreement, this discrepancy could be easily corrected by 
adjusting the efficiency of H 2 FD. 


A likely explanation for these differences is generally lower CO and CO 2 ice abundances 
toward high-mass objects that experience rapid contraction and darkening. Observations 
of protostars indicate that this indeed may be the case, with high-mass objects having 


roughly a half of t 


le carbon oxide inv entory observed in low- or intermediate-mass objects 


Gibb et ah 

2004; 

Oberg et ah 

2011bj 


2011bl ). The abundances of solid H20:C0 and H20:C02 


ratios for Models L and M approximately agree to these observations. Thus, it can be said 
that the presented conclusions on selective desorption mechanisms are valid for a wider 
range of prestellar objects, not just low-mass cores. 


3.5.1. The chemistry of COMs 

As stated in section 13.11 the model produces interesting results regarding complex 
organic molecules. COMs are efficiently produced in the sublayers by the photoprocess, 
and are released into the gas by photodesorption of species in shallow ice layers directly 
below the surface fsection 12.31) . Hydrogen-poor molecules, such as formic acid, formamide 
HCONH 2 , formaldehyde CH 2 O, glyoxal (CHO) 2 , and other, heavier COMs may contain 
a significant part of oxygen and carbon reservoir. For the complete model, this produces 
a relatively high gas-phase abundances of many organic species between approximately 
0.73 and 0.96Myr, or Ay 2.4 to lOmag. The first number indicates the formation of first 
sublayers and the onset of subsurface PDI chemistry. The second number signals the end 
of the ISRF photons as a major source of desorption and dissociation, and the initiation 
of the rapid accretion phase, which prevents desorption of species whose formation time 
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Fig. 8.— Gas-phase (top panel) and ice (bottom panel) abundances, relative to hydrogen, 
of the complete Model L for organic species observed in quiescent dark clouds. 
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Fig. 9.— Gas-phase abundances, relative to hydrogen, of selected mantle-produced species 
for variants of Model K with a = 0.01 (top panel) and a = 0.06 (bottom panel). 
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is longer that the time required to ‘bury’ them in the mantle. Additionally, during the 
rapid freeze-out, thickness of the surface layer exceeds 2ML, at which point desorption 
from the sublayers is prohibited (section [53]) • As stated in section IXTl the model produces 
interesting results regarding complex organic molecules. COMs are efficiently produced 
in the sublayers by the photoprocess, and are released into the gas by photodesorption 
of species in shallow ice layers directly below the surface (section 12.3p . Hydrogen-poor 
molecules, such as formic acid, formamide HCONH 2 , formaldehyde CH 2 O, glyoxal (CHO) 2 , 
and other, heavier COMs may contain a signihcant part of oxygen and carbon reservoir. For 
the complete model, this produces a relatively high gas-phase abundances of many organic 
species between approximately 0.73 and 0.96Myr, or Ay 2.4 to lOmag. The first number 
indicates the formation of first sublayers and the onset of subsurface PDI chemistry. The 
second number signals the end of the ISRF photons as a major source of desorption and 
dissociation, and the initiation of the rapid accretion phase, which prevents desorption of 
species whose formation time is longer that the time required to ‘bury’ them in the mantle. 
Additionally, during the rapid freeze-out, thickness of the surface layer exceeds 2ML, at 
which point desorption from the sublayers is prohibited 1 section l2.3p . 


The reac t ion ne twor k has already been described in detail by 


Garrod fc Herbst 


Garrod et al. 

(2008 

) and 

Laas et al. 

(2011) 


( 120061 ); 


(1201 ll ). These authors inv estigate the formation of 


COMs during the gas warm-up phase in the protostellar stage. 


Vasvunin &: Herbst 


fl2013bl) 


investigate the gas-grain chemistry of COMs in cold cores. They use reactive desorption 
as the main means for heavy molecule ejection into the gas phase. They consider a clump 


dense gas (lO^cm existing in steady physical conditions 


years. The best-fit results obtained by fjVasvunin fc Herbst 


or several hundred thousand 


2013bl) are for model with a 


very efficient reactive desorption process - 10% of reaction products go to the gas phase, 
and all ice species are affected by surface reactions (no bulk ice). Because of this, reactive 
desorption is attributed to all ice species and is summarily more effective roughly by a factor 
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of 10^ than in Models J to M. 


Vasvunin &: Herbst 


fl2013bl) conclude that the assistance of 


other processes is required to reproduce gaseous COM abundances in protostellar cores. 


The research presented here has shown that reactive desorption is very ineffective in 
the case of COMs. Meanwhile, surface synthesis of COMs is made more productive, because 
the large molecules remain on the grains. Moreover, it has been shown that photoprocessing 
of thin ices at low Ay, followed by photodesorption from shallow subsurface layers, can be a 
plausible source for COMs in the gas phase. This process is eff icient at much lower densiti es 
and shorter evolutionary time-scales than those considered by 


Vasvnnin fc Herbst 


(12013^. 


Figure [ 8 ] shows the calculated Model L abundances, relative to H, for a selection species 
that have been observed in dark clouds. These species include methanol CH 3 OH, formic 
acid HCOOH, propynal HC 2 CHO, acetaldehyde CH 3 CHO, methyl formate HCOOCH 3 , 
ketene CH 2 CO, dimethyl ether CH 3 OCH 3 , methoxy (CH 3 O), and formyl (HCO) radicals 
(see Table [ 8 ]). All the observed abundances fall in the vicinity of 10“^° — 10“^^, and the 
Standard model produced abundances in or near this range. The changes in desorption 
efficiencies in the complete model result in a poorer agreement with observations for several 
of the species. 


In the following paragraphs a short analysis on the Model L chemistry for significant 
organic compounds, observed in quiescent cores, is presented. Table [ 8 ] shows a comparison of 
observational and calculated abundances of organic species in dark clouds. The conditions 
and history of the regions sampled by observations are unknown. An agreement with 
observations can be claimed in cases, where the calculated maximum abundance equals or 
exceeds the observed value. 


There is a dip for the ice abundance curves of radical species HCO and CH 3 O, 
associated with the weakening of the ISRF. A similar behavior can be observed for most 
hydrogenated radical species. All radicals then experience a steep upward trend in the hnal 
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stages, associated with the overall freeze-out of molecules. 

Formaldehyde. H 2 CO is synthesized in ice in signihcant amounts (0.1% of H 2 O at 
IMyr) via barrierless double-hydrogenation of CO by atomic H, although other pathways 
exist. It can be photodissociated back to CO or HCO, or converted to HCO by its reaction 
with the abundant OH radical. Other surface binary reactions of H 2 CO have relatively high 
energy barriers and are ineffective. 

Formic acid. A barrierless reaction HCO-I-OH is the main pathway for the formation of 
HCOOH. Thus, formic acid is the end product of a barrierless reaction chain that involves 
CO and the abundant radicals H and OH. As such, HCOOH is produced in vast amounts in 
the sublayers, where OH is readily available. In Model L, the abundance of formic acid in ice 
reaches 1% relative to that of H 2 O. This is a significantly better result than the unrealistic 
15% for the Standard model. HCOOH is overproduced in the Standard model because of 
the prolonged PDI period, when CO is converted into HCOOH and other compounds. In 
turn, this period is so long because of the inefficient desorption in that model. 

Other species formed in the mantle. Methyl formate HCOOCH 3 has a low ice 
abundance of 6 x 10“^ relative to H 2 O at IMyr. Its production requires the protonation 
of formaldehyde - a reaction with a barrier. The same holds true for dimethyl ether 
CH 3 OCH 3 , whose synthesis does not involve any major radical species and thus is even 
lower than that of other COMs. Acetaldehyde CH 3 CHO is able to form via the CH 3 -I-HCO 
reaction. The synthesis of these two radicals involves no barriers (similarly to methanol), 
and CH 3 CHO reaches somewhat higher abundance than other complex species. 

Surface species. Several of the species in Table [8] can be associated exclusively with 
surface or gas-surface processes. The formation of CH 3 OH occurs via an interaction of 
gas and surface chemistry. The hydrogenation of formaldehyde is hampered by activation 
barriers, and most of CH 3 OH forms via the CH 3 -I-OH reaction. This is a barrierless pathway 
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that requires the successive hydrogenation of C with three H atoms. The radicals C, CH, 
CH 2 , and CH 3 can also be available from the accretion of gas-phase species. The resulting 
production rate in ice for CH 3 OH is low. It can be signihc antly higher in protoste llar 


envelopes, where reactions with barriers become important (IGarrod fc Herbst 


200611 . The 


metoxy radical CH 3 O almost exclusively originates from the photodissociation of methanol. 


Ketene H 2 CCO and propynal HCCCHO are also associated with surface chemistry. 
Their synthesis starts with radicals that accrete directly from the gas, e.g., H, C, O, and C 2 . 
Although barrierless, the synthesis involves a number of steps, and has a limited efficiency. 
The abovementioned radicals are not abundant in the sublayers. Because of this, H 2 CCO 
and HCCCHO do not reach particularly high abundances in ice. However, the gas-phase 
abundance of these compounds is relatively high because of their concentration on the 
surface. This is especially obvious for the early core evolution stages (hgurelHl). 

Formic acid and formamide are produced in PDIs in vast amounts. These species 
form from radicals that are readily available in subsurface ice - OH, HCO, and NH 2 . The 
overproduction largely arises from the fact that most current reaction networks, including 
the one employed here, are designed for surface chemistry. Reactions for organic species 
with radicals that are almost non-existent on surface but abundant in subsurface layers - 
NH, NH 2 , OH, HCO, and others - have largely been omitted in the network. These radicals 
arise from the direct photodissociation of the most common ice constituents, and should 
have a tremendous importance in the synth esis of complex molecules . A similar conclusion 


on subsurface ice chemistry was reached in 

Kalvans & Shmeld 


201 

G). The importance of 

such radicals in mantle reactions has been recognized by 

Garroc 


>013a 

), who included 


hydrogen abstraction reactions by OH in his three-phase model. 


CarrodI fj2ni3al ) successfully modeled the abundances of organic molecules during the 
warm-up and eventual evaporation of the (circumstellar) ices in the envelope of a protostar. 
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In such conditions, cosmic-ray-induced photons are found to be important for radical 
generation in ices, instead of ISRF photons. 

The main result from the Standard model remains valid also in the complete model - 
organic molecules can be produced in PDIs in signihcant amounts, and their non-thermal 
release into the gas can be the responsible for the gaseous presence of at least several 
organic compounds in dark clouds. The two other complete models - J and K - produce 
gas-phase COM abundances within the same order of magnitude. 

Because of its peculiar selectivity, direct reactive desorption may alter COM abundances 
in gas and ice. This mechanism has a dual effect - more efficient desorption even more 
delays the formation of subsurface ice and the synthesis of COMs in the sublayers. For 
example, if the value of a is taken 0.1 instead of 0.03 for Model L, the first sublayer forms 
only at 0.855Myr instead of 0.730Myr. However, small molecules are much more effectively 
desorbed, and the overall proportion of COMs in ice is higher, once the first sublayer has 
been initiated. 

Variations of Models J, K, and L with a values in the range of 0.01-0.06 result in 
abundance differences up to several tens of per cent for the organic species. Figure [9] shows 
an example of gas-phase abundances with two different a values for selected species that are 
produced mostly in subsurface ice. The behavior of individual species can be highly specific. 
A general trend is a slight increase for gas-phase abundances of the heavier COMs at lower 
a values. Reactive desorption is highly dependent on the number of reactions involved in 
the formation of each species, and on the value of the parameter /^d for each particular 
reaction. This is unlike the derivations (Models B-D) of the Standard model, where a 
higher a meant higher abundances of COMs. For the complete model, other desorption 
mechanisms hamper the formation of the sublayers, and any increase in the efficiency of 
reactive desorption delay the formation of subsurface ice even more, thus delaying the 
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photo-synthesis of COMs. The overall result is that, counter-intuitively, more effective 
reactive desorption produces lower gas-phase abundances of COMs and other subsurface 
species. This also holds true for the solid phase. 

Because complex species have similar assumed Tph, photodesorption has no significant 
effect on COM abundances. Derivations of the Standard model have negligible differences 
on the order of a few per cent at most (section l3.2p . The effect of photodesorption in the 
complete models is even lower because sublayer ices are exposed to the ISRF for a short 
time, only. 


3.5.2. Other ice components 


Figure [To] shows the abundance of selected species in the ice and the gas phase. Ice 
species can be broadly divided into two categories - those produced on the surface (H 2 O, 
CH 3 OH, CO 2 , H 2 CO, H 2 S), and those produced in the ice mantle (O 2 , HCOOH, SO, SO 2 , 
OCS, H 2 O 2 ). During the late stages, a third category appears - the molecules accreted 
directly from the gas phase (CO, N 2 , HNC). Naturally, intermediate sub-categories exist 
with signihcant contributions from two or three categories (HCN, O 2 ). 


Sulfur chemistry is dominated by surface production of H 2 S. SO is the second most 
important sulfur molecule in ice with a final abundance of 42% relative to that of H 2 S. 


Oxides and OCS appear during late stages in the subla yers and reach a 


Dundances of 0.2 


and 0 . 02 % rela 

dve 

0 H 9 S, respectively. ] 

(Booeert et ah 

1997; 

Zasowski et ah 

0 

0 

CM 


1997 1 and SO 2 


2009!) are the only solid sulfur compounds observed 


towards young stellar objects, a possible implicati on is that subsurface che mistry plays a 


significant role in grain surface processes (see also 


Kalvans fc Shmeld 


201011 . 


The calculation results show an overabundance of solid ammonia, which is typically 
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Fig. 10.— Gas-phase (top panel) and ice (bottom panel) abundances, relative to hydrogen, 
of the complete Model L for selected species. 
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not observed in dark clouds fjOberg et alJl2011b|) . The H 20 :NH 3 ice abundance ratio is 


17% for Model 


and 20 % for t 


similar results flGarrod fc Paulv 


re Standard model. Relevant mod els by other authors yield 

2013a). The formation of NH 3 


20111 : 


Vasvunin fc Herbst 


occurs simultaneously with H 2 O, and via similar reactions. Because these two species have 
very similar Ed, they cannot be separated by selective desorption. A probable cause for 
such results is an underproduction of gas-phase N 2 . This problem is only partially solved 
by the inclusion of N 2 shielding. A dedicated investigation of nitrogen gas-grain chemistry 
may be required to hnd a solution for this problem. Ammoni um ion NHt can be a m ajor 


reservoir for nitrogen in ice, as suggested in previous studies fjSchutte fc Khanna 


200811 . 


Hydrogen peroxide has be en observed in the int erstellar medium with a high abundance 


of 10 relat ive to hydrogen 


chemistry by 


Du et ah 


(iBergman et ah 


201 ll) . This has been attributed to surface 


(1201211 . with reactive desorption as the main mechanism. These 


authors also investigate in detail the formation mechanism of H 2 O 2 . Th is model was then 


verihed by the observation of the O 2 H radical with a similar abundance fjParise et al 


2 OI 2 I). 


In the complete model presented here, gaseous H 2 O 2 is produced with relative 
abundances up to 10“®, in agreement with observations. Figure 13.5.21 shows a double 
production peak. H 2 O 2 is produced by surface reactions, while subsequent hydrogenation 
gradually lowers its abundance. A second peak is observed with the onset of active mantle 
photochemistry. The model fails to produce gas-phase O 2 H at rela tive abundances higher 
than 10“^^. Differences between this model an the one employed by 


Du et al. 


(1201211 include 


the activation energy for the reaction H -f O 2 (1200 and 600K, respectively) and different 
physical conditions. These likely do not produce signihcant changes for O 2 H gas-phase 
abundance. 


by 


Most importa n tly, reactive desorption is m ore effective by a factor of 10^ in the models 


Du et al. 


(I2QI2I); 


Vasvunin fc Herbst 


(l2013bfl . This is because of their high a = 0.1 value. 
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and because the hundred-monolayer-thick bulk mantle was not considered by these authors. 


The effectiveness of reactive desorption has to be ca refully eva l uatec 


two-phase models. It is possible that, if the model by 


Du et ah 


before application in 


(120121) considered subsurface 


mantle, then it could be difficult for this model to attain high gas abundances for species 
produced on surface. 


The current paper 
abundance observed by 


certainly cannot of fer a plausible explanation for the high O 2 H 


Parise et al 


fl2012h . The molecular cloud core in consideration, p 


Oph A is notable with observations o f molecular oxygen and probably is in a transient 


evolutionary stage flLiseau et al 


20121). A perturbed state is probably supported by the 


relatively high gas temperature of 20-30K, uncharacte ristic for c 


of several protostars in the direct vicinity of the core fiDu et al. 


(Du et al. 

2012 ; 

Liseau et al. 

2012 ) 


3.5.3. The composition of the sublayers 


In a fully formed ice at IMyr, sublayer 3 has a Il20:C0:C02:NIl3:N2 ratio of 
60:15:10:12:4, respectively. For sublayer 2 this ratio is 39:41:7:6:8, for sublayer 1 - 
22:54:7:2:13, and for the surface layer - 7:50:12:1:21. The outer sublayer 1 is thinner - 
«i40ML instead of 60 for a ‘full’ sublayer. Figure ITT] shows the evolution of water and 
carbon oxide abundances in the sublayers. Sublayers 1 and 2 represent the t wo observed 


modes of solid CO - in polar and in apolar ice matrix fjSandford et al 


199l|). 


1988 : 


Tielens et al 


Bergin et al. 


( 2005 ) conclude that towards the field star Elias 16 {Ay = 24.1) 15% of 


CO 2 resides in apolar ice. For Model L at the corresponding Ay the total ice thickness 
is 92.6ML, and only sublayer 3 is ‘full’ with 60.6 MLs of ice in it. 74% of CO 2 reside in 
sublayer 3, which is water-dominated (Il20:CO 100:25), and the remainder in other layers. 
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Fig. 11.— Model L results: abundances of water and carbon oxides in monolayers per 
sublayer, relative to total H, during the last lOOkyr of cloud contraction. The curved pattern 
for surface (upper panel) abundances is an artifact because of the surface-sublayer transition 
as described in section 12.4.21 Note the different abscissa scale for the surface layer. The 
composition of sublayers also reflects the evolution of surface composition over time. 
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where H20:C0 is roughly 1;1. These results probably are in an approximate agreement 
with observations. The composition of the surface layer indicates that the apolar CO most 
likely resides in the topmost layers if interstellar ices. 


4. Conclusions 


A three-phase model has been presented with several features that may aid in a better 
understanding of interstellar ice-related chemistry and physics. The concept of sublayers 
facilitates the modeling of chemical reactions in ice. This approach also has the potential 
to simulate protostellar ice evaporation in a com plex way, as the sublayers are gradually 


exposed to the surface (see 


Favolle et ah 


2nilal ). This will be used in subsequent papers. 


The transition of molecules from the surface to the mantle with a hnite rate allows the 
possibility to simulate the compaction of an initially porous ice (although this has little 
effect on ice composition). To a limited extent, desorption from the subsurface layers is 
permitted, which may help to explain the observed abundances of at least some organic 
molecules in dark clouds. Finally, the research suggests the existence of ‘photon-dominated 
ices’ in weakly shielded regions of molecular clouds. Basically, shallow ice layers, isolated 
from the surface, may experience intensive processing by interstellar photons. This gives 
rise to a peculiar chemistry on grains during early stages of core contraction - or in the 
exterior part of the cores. 


The influence of three desorption mechanisms on ice composition (the H20:C0:C02 
ratio) has been investigated in detail. Table [9] shows the summary of the specific models 
K and L are the ‘hnal’ models that yield the best agreement with observations for the 
Ay-dependent H20:C0:C02 ratio 1 section 1^3]) . 

Photodesorption is fairly well constrained by laboratory results (Table [3]) and the 
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adoption of Fph values for the ISRF photons that are in agreement with experiments is 
essential for a proper reproduction of ice accumulation in early core contraction stages. 
The values suggested here are 3 x 10“^ for ISRF photons and 2 x 10“^ for cosmic-ray 


i nduced photons, instead o f a single value o : 


fIVasynnin fc Herbst 


2013a 


bl; 


Garrod 


2013a : 


10 used in many current researc. 


Gerner et ah 


2014 : 


1 pap ers 


Ghang fc Herbst 


20141 1 


For reactive desorption, it 


las sometimes been assumed that all binary r eactions 


have similar a) efficiencies flGarrod et ah 


20071 : 


Vasvunin fc Herbst 


2013b 1. w: 


other authors use rea ction- and molecule-specific desorption efficiencies flDu et ah 


rile 


2012 


Reboussin et ah! 


2014J) . The latter approach is the one used in the present study. No 


detailed analysis of selective reactive desorption could be found in the literature. The 
study of this mechanism in section 13.31 has yielded some unexpected results. Reactive 
desorption is inefficient for atom exchange reactions, characteristic for the formation of 
large molecules. These reactions typically do not produce much heat. Instead, the highly 
exothermic atom-addition reactions involved in the formation of simple species result in 
effective desorption for water, ammonia, and carbon dioxide. 


Reactive desorption has a diverse effect on the abundances of COMs. Efficient 
desorption (higher a) results in the accumulation of atoms in heavy molecules, because 
light species accumulate slower. This results in a thinner ice during early stages, and 
first sublayers appear later, which may signihcantly shorten the time available for organic 
synthesis in the sublayers. The latter effect is more important, because the effectiveness 
of PDI chemistry is strongly bound to ice thickness at low Ay values. In any case, a 
reaction-specific approach on reactive desorption is essential for modeling of COM synthesis 
on interstellar dust grains. 


The third mechanism in consideration was indirect reactive desorption, with desorption 
arising from the H-f H reaction as its most important representation. It has an efficiency 
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that is very poorly constrained by experiments or theory. It was fonnd that this mechanism 
is probably responsible for the observed H20:C0:C02 ratio in interstellar ices. As a 
requirement for this, H 2 FD has to be efficient for CO and CO 2 , and negligible for water. 
Currently, as far as I know, there are no data that can either confirm or deny such 
a conclusion. The estimated efficiency parameter /nafd likely has to be in the range 
10“® — 10“^ desorbed molecules per accreted H atom for the CO molecule in order to 
produce a good agreement (±50%) with observations. These /nafd values do not counter the 
existing knowledge in the sense that they are two to four orders of magnitude lower that 
the efficiency of direct reactive desorption. Naturally, H 2 FD is bound to the abundance of 
H atoms in the cloud, which might be a cause of uncertainty for fn^fd- 

There are three general conclusions on ice accumulation. First, the onset of ice 
accumulation onto grains is mostly governed by desorption by interstellar photons, 
although direct and indirect reactive desorption also have an effect. Second, the observed 
Ay-dependent H20:C0:C02 ratio can be reproduced by a mechanism that is sufficiently 
effective and has a sharp difference in its yield for carbon oxides and water. Desorption by 
the H±H reaction on grains is a likely candidate. Third, efficient desorption during early 
stages of cloud contraction result in nominally slightly thicker ices in the dense cloud core. 
This is because the longer an ice layer is exposed to the ISRF, more multi-atom species are 
synthesized, and the total number of molecules is lower. 

The three sublayers of subsurface ice represent three chemically distinct ice components 
at a stage when the freeze-out has ended. Starting from ice layer near the grain core itself, 
sublayer 3 consists of H 2 O with an admixture of NH3 and CO 2 , sublayer 2 is a 1:1 mixture 
of H 2 O and CO, and sublayer 1 is a 2:5 H20:C0 mixture with an addition of N 2 and CO 2 . 
The surface is largely covered with CO. These results are in compliance with observations 
of polar (water excess) and apolar CO and CO 2 ices fsection r3.5.3p . 
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As noted in section 13.4.21 all the models considered in this paper show an apparent 
depletion of CO ice at its respective Ath when compared to H 2 O and CO 2 . This result 
pertains also for the complete models (Table [7]). Because the A^^-dependent H20:C0:C02 
ratio for Models K and L matches the observations, and because H 2 O and CO 2 abundances 
at their respective Ath are quite similar, it can be possible to deduce a new detection 
threshold Ay for solid CO. Based on these results, an estimate of Ath for CO is 10.5 from 
Model K and S.lmag from Model L. Only the latter value li es within the error margin of 


the initially estimated Ath = 6.8 ± 1.6 from observations by 


Bergin et al. 


( 2 OO 5 I). Under 


the assumption that ice species have similar abundances at their respective Ath, the model 
results suggest that the threshold Ay for CO ice probably is in the range of 8-10.5mag, 
with the lower value being supported by currently published observations. 


The higher spread for the calculated-to-observed CO ice abundance ratio (Table [7]) 
probably arises because, unlike water and CO 2 , CO basically does not accumulate in ice 
until very late times. This means that a huge mass of CO accretes onto the grains on a very 
short time-scale, when the cloud has become sufficiently dark and dense. Consequently, 
relatively minor physical perturbations during late evolutionary stages of the core, caused 
by, e.g., nearby stars or outflows may speed-up or delay the CO, N 2 , and O 2 accretion peak, 
producing the observed spread in abundances relative to H 2 O, CO 2 , and NH3. The latter 
molecules can be expected to be less prone to such temporal perturbations because they 
begin to accumulate much earlier. 
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Table 7: Comparison of calculation results the complete model, variants J, K, L, and M (see 
text and Table [3 for an explanation). 




J 


K 


L 


M 

Av 

co^ 

C02^ 

co^ 

C02^ 

co^ 

C02^ 

co^ 

C02^ 

6.3 


1.0 


1.3 


1.2 


0.9 

10 

1.0 

0.9 

0.5 

1.1 

0.8 

1.0 

0.4 

0.8 

11.7 

2.0 

0.9 

1.1 

1.2 

1.6 

1.0 

0.8 

0.8 

15.3 

1.1 

0.9 

0.7 

1.0 

0.9 

0.9 

0.4 

0.6 

17.8 

0.8 

0.8 

0.5 

0.9 

0.6 

0.8 

0.3 

0.5 

20.9 

1.2 

0.9 

0.8 

1.0 

1.0 

0.9 

0.4 

0.5 

21.5 

1.7 


1.2 


1.5 


0.6 


24.1 

1.8 

0.7 

1.3 

0.8 

1.6 

0.7 

0.6 

0.4 


H20:C0:C02 (ML) 

H20:C0:C02 (ML) 

H20:C0:C02 (ML) 

H20:C0:C02 (ML) 

Ath 

3.5 : 

1.9 : 2.4 

3.6 : 

1.0 : 3.8 

2.5 : 

1.3 : 2.2 

3.7 : 

0.4 : 1.6 
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Table 8: Comparison of observed and calculated (maximum) abundances of complex organic 
molecules in quiescent cloud cores. 


Species 

Observations, cm 

Bl-b L1689b 

3 

Model L, cm ^ 

Gas^ Ice^ 

H 2 CO 

4.0E-10 

M3 1.3E-09 


3.4E-10 

1.2E-07 

HCO 

1.8E-11 

C® 



2.4E-10 

2.2E-09 

HCOOH 

l.OE-11 

C 



4.2E-08 

1.4E-06 

HCOOCH 3 

2.0E-11 

C 

7.4E-10 

B 

9.3E-13 

7.7E-11 

CH 3 OH 

3.1E-09 

06 



1.4E-10 

3.9E-09 

CH 3 O 

4.7E-12 

c 



2.4E-11 

1.6E-12 

H 2 CCO 

1.3E-11 

c 

2.0E-10 

B 

7.8E-12 

9.2E-13 

CH 3 CHO 

l.OE-11 

c 

1.7E-10 

B 

2.4E-13 

1.3E-10 

HCCCHO 

3.6E-12 

c 



8.2E-12 

2.5E-12 


^Maximum abundance 


^At t = IMyr 


" Marcelino et al. 

2005 


Bacmann et al. 

2012 ) 

' Cernicharo et al 

( 

201 

1) 

^ Obere et al. (2010 

) 
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Table 9: List of desorption parameters employed in models. 



Pliotodes. 

a, reactive 

HsFD 

H 2 FD Ed 

Model 

yield 

desorption 

efficiency 

threshold (K) 

Standard 

0.001 

N/A^ 

0 

N/A 

A 

uniform^ 

St.3 

St. 

N/A 

B 

St. 

0.01 

St. 

N/A 

C 

St. 

0.03 

St. 

N/A 

D 

St. 

0.06 

St. 

N/A 

E 

St. 

St. 

4 X 10"® 

1210 

F 

St. 

St. 

2 X 10-® 

2600 

G 

St. 

St. 

2 X 10-® 

6000 

H 

St. 

St. 

RRK 

6000 

I 

St. 

St. 

emp.^ 

none 

J 

uniform 

0.03 

4 X 10-® 

2600 

K 

uniform 

0.03 

RRK 

6000 

L 

uniform 

0.03 

emp. 

none 

M® 

uniform 

0.03 

RRK 

6000 


^Reactive desorption efficiency assumed 1% 
^0.002 for secondary and 0.003 for ISRF photons 
^As in Standard model 
^empiric dependent relation 
®Short (0.2Myr) core contraction time 
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